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METHOD FOR CALCULATING THE AERODYNAMIC LOADING ON AN OSCILLATING FINITE 

WING IN SUBSONIC AND SONIC FLOW 1 

By Habby L. Runyan and Donald S. Woolston 


SUMMARY 

This report presents a method for determining the air forces 
on an oscillating finite wing of general plan form in subsonic 
flow including the limiting case of sonic flow. The method 
utilizes some of the concepts developed by Falkner ( British 
R. & M. No. 1910) for steady lifting-surface theory. The load- 
ing on the wing is assumed to be given by a series containing 
unknown coefficients which satisfies various boundary con- 
ditions at the edges. The reguired integrations are performed 
by approximate means, and a set of simultaneous equations in 
terms of the coefficients in the loading series is obtained. 
Solution of this set of equations then gives the loading coefficients. 
The method is applied to rectangular and delta wings and com- 
parison is made with existing theory. A sample calculation 
is given in an appendix. 

INTRODUCTION 

The analytical determination of the air forces on oscillat- 
ing wings has been a continuing problem for oyer 30 years. 
Most of the effort has been directed toward the determina- 
tion of the forces on wings in two-dimensional flow, and re- 
sults have been obtained for a complete range of Mach 
numbers, both subsonic and supersonic. For finite wings, 
however, the analytical work is still in a state of develop- 
ment for all speed ranges. The main effort has been di- 
rected toward the incompressible case (for example, refs. 1 
to 8) and the supersonic field (refs. 9 to 14) leaving an im- 
portant area virtually untouched, namely, the subsonic range 
of Mach numbers between 0 and 1. Among the few studies 
that have been made in this range of Mach numbers have 
been those of Merbt and Landahl (ref. 15), who extended 
some of the work at very low aspect ratio to include the ef- 
fect of Mach number, and the work of Yoss, Zartarian, and 
Hsu (ref. 16). In addition, Reissner (ref. 17) made a pre- 
liminary assessment of the use of Mathieu functions for 
subsonic compressible flow, and W. P. Jones (ref. 18) and 
Lehrian (ref. 19) studied the cases for larger aspect ratios. 
The present report deals with some recent efforts which 
have been made toward filling this gap, namely the develop- 
ment of a wing-surface theory for subsonic and sonic flow. 

In the stationary two-dimensional case, the effects of com- 
pressibility may readily be obtained from results of incom- 
pressible theory by application of simple transformation or 
correction factors such as the well-known Prandtl-Glauert 
factor. In the stationary three-dimensional case, the effect 


of compressibility may be treated as an extension of the two- 
dimensional procedure in a simple maimer as proposed by 
Gothert (ref. 20). In a further extension for the oscillating 
compressible case, Miles (ref. 21) has proposed a method in- 
volving several transformation steps for correcting the re- 
sults of incompressible theory to include the effects of com- 
pressibility for a finite wing. The method is, however, re- 
stricted to low values of the reduced frequency. It there- 
fore appears to be necessary to deal with the boundary-value 
problem directly in order to provide a method of general 
applicability which may he used, for example, in flutter work. 

The linearized boundary-value problem for the oscillating 
finite wing in a compressible flow may be approached from 
two points of view. One approach involves the transforma- 
tion of the governing differential equation by a suitable 
choice of coordinates followed by the use of the classical 
method of separation of variables. Solutions are then found 
in terms of series of orthogonal functions. This method was 
used by Schade and Krienes (ref. 22) in obtaining solutions 
for the oscillating circular plate and has recently been gen- 
eralized by Kuessner (ref. 23). It has been pointed out that 
there is a definite limit to the number of generalized coordi- 
nate systems appropriate to the wave equation and to the 
corresponding orthogonal functions available. The method 
is, therefore, limited to specialized plan forms such as the 
circular or the elliptic plate. 

The second approach involves a direct consideration of 
the integral equation relating pressure and downwash. The 
integral equation may be derived from the standpoint of 
either the velocity potential or the acceleration potential. 
Inasmuch as one is primarily interested in the pressures on the 
wing surface, it seems more direct to use the form of the 
integral equation associated with the pressure, namely, the 
acceleration potential approach; therefore, this procedure is 
used in the present report. 

Basic to the solution of the integral equation is the eval- 
uation of its kernel, a function which represents the down- 
wash at a point in the plane of the wing due to a unit loading. 
Only recently (ref. 24), the kernel function for oscillating 
finite wings in subsonic and sonic flow has been reduced to a 
form which can he conveniently evaluated. This back- 
ground of available information on the kernel has led to the 
consideration of numerical approaches for handling the 
integral equation in the development of general lifting- 
surface methods which could be used for any speed range. 


> Supersedes NAOA Teohnlcal Note 3631 by Harry L. Banyan and Donald S. Woolston, 1858. 
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A similar integral equation appears in the steady case, but 
involves a much simpler kernel function and a number of 
approximate solutions have been successfully developed. 
In selecting an appro xim ate method for handling the integral 
equation in the oscillatory case, the various procedures 
developed for the steady case were examined. In con- 
sidering these procedures, it was kept in mind that a lifting- 
surface theory which is to be used for flutter calculations 
must meet certain requirements beyond those usually con- 
sidered in the steady case. One essential requirement is that 
not only the lift but also the moment must be accurately 
predicted. In several of the methods now being used for 
the steady case, the wing is essentially replaced by a single 
line at the X-chord position. This placement essentially 
fixes the center of pressure at the X chord, and thus the 
moment is constrained to become a function only of the lift 
force. In the case of an airplane of normal configuration 
with a tail placed a considerable distance from the wing, the 
contribution of the wing pitching moment is very small 
compared with the pitching moment due to the tail. How- 
ever, for tailless configurations and for aeroelastic problems 
such as flutter, the moment on the wing is of predominant 
importance and must be accurately predicted. 

In addition to the previously mentioned "requirement of 
predicting the moment, the method should be easily adapt- 
able to the calculation of the loading on oscillating flexible 
wings and, in addition, should take into account the effects 
of compressibility. It appeared that, among others, a 
method developed for the steady incompressible case by 
Falkner (ref. 25) could be extended to include these require- 
ments and this method was selected as the basis for the 
present investigation. 

Among other methods considered was a multiple-line 
approach suggested by the procedure of Schlichting and 
Kahlert (ref. 26) for the case of steady flow. The procedure 
is thought to be of interest in that it indicates a simpler 
alternative -«neans of handling the lifting-surface problems. 
An extension of this method to the oscillatory case is also 
made in the present investigation. 

The primary aim of the present investigation has been the 
development of a lifting-surface method for calculating the 
forces on a wing of any plan form which is harmonically 
oscillating in a subsonic or sonic flow. The method uses 
some of the concepts of Falkner in an extension to the com- 
pressible oscillating case. First, a brief description of the 
basic integral equation and some of the possible means of 
solution are given. Then, the basic theory of the method is 
given in detail in the body of the report, and the results of 
some calculations are shown for rectangular and delta 
wings. In appendix A, the numerical details of calculating 
the forces on a rectangular wing are given as an example. In 
appendix B, a treatment is given for certain integrals which 
arise in the method and which contain singularities. In 
appendix C the expression for the pitching moment for a 
delta wing is obtained. Finally, in appendix D, a descrip- 
tion and some results of applying a multiple-line approach, 
based on a procedure used by Schlichting and Kahlert for 
steady flow, are given. 


SYMBOLS 


A 

aspect ratio 

a 

axis of rotation measured from 


midchord in terms of half-chord) 
positive aft 

a>im 

coefficients in expression for loading 

b 

half-chord, ft 

&!,<, 

lift coefficient associated with a 


degree of freedom g 

Cm, 5 

moment coefficient associated with 


a degree of freedom g 

c 

chord, ft 

s 

II 

1 

!PI*> 


F'tr 

nondimensional downwash factor 


defined by equation (29) 

/*(£) 

chordwise loading function 



Of o, Ofn 
Qmiyi) >9m (/*) 

l)] 

chordwise replacement loads 
spanwise loading functions 

kernel of three-dimensional in- 

tegral equation 

K(M,Z) 

kernel of two-dimensional integral 


equation, 

~E{M,Z)= 


k[.M, yix-S^iy-viJdv 

s 

nondimensional three-dimensional 


kernel defined by equation (28) 

k 

reduced-frequency parameter 


based on half-chord, bu/V 

k 

reduced-frequency parameter 


based on «, 

La . 

total lift on whig, associated with 


a degree of freedom g, lb 

1 

reference length 

M 

Mach number 

M, 

moment about a, associated with 


a degree of freedom g, ft-lb 

n,m,j 

integers 

2 

dynamic pressure, ^ pV 1 

#i=V(z-£) s +/S 2 (y+ D 1 

i? 2 =y(x-?) s +is 2 (r-i ) 2 

S 

area of wing surface 

8 

semispan of wing, ft 

t 

time 

V 

w(x,y,f) 

stream velocity 

vertical induced velocity or down- 

wash 

x,y,z,t,v 

rectangular coordinates 

Xo= (x—x,) 

control-point location 

Xr 

chordwise coordinate for replace- 


ment loads 
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Vs 


Vo=y-v 

y'L.y'n 


Z—k(;Xg—£) or k(x t —x,) 
a 

/s=vr 

r a,u) 

A P,(m) 

a 

s 


X 

p 

4- 

<t> 

(J> 

Subscript: 

2 


distance from control point to 
center of segment 

distances from a control point to 
midpoints of load segments lo- 
cated, respectively, to the left 
and right of the wing midspan, 
referred to e 

generalized' coordinate of an ang- 
ular displacement 

vorticity distribution 
loading function defined by equa- 
tion (26) 

pressure difference 
thickness ratio 

semispan of spanwise load seg- 
ment, s/20, ft 

angular ehordwise coordinate, 
x=— cos 6 

variable of integration 
air density 

angular spanwise coordinate, 

1 J = — cos <t> 

phase angle between lift or mo- 
ment and position 
circular frequency, radians/sec 

denotes a degree of freedom 
(for pitch, q—a) 


DISCUSSION OF INTEGRAL EQUATION AND METHODS OF 

SOLUTION 


The determination of the aerodynamic forces on an 
oscillating wing is a problem which, in general, must be 
handled by approximate or iterative procedures. In the 
two-dimensional case one approach to the problem has been 
the transformation of the governing differential equation 
and its associated boundary condition^ to an integral equa- 
tion for which approximate solutions can be found. It is 
interesting to note that for this two-dimensional case, an 
iterative solution based on the integral-equation approach 
preceded a more exact solution based on a series of Mathieu 
functions by about ten years. 

In the case of wings of finite span, a transformation 
similar to that made in the two-dimensional case and leading 
to an integral equation can be made. The integral equation 
may be derived on the basis of the velocity potential or on 
the basis of the acceleration potential. The advantage of 
the use of the acceleration potential rests in the fact that 
the trailing surface of discontinuity is not explicitly intro- 
duced into the problem and the integral equation relates 
directly the downwash and the pressure difference on the 
wing. The disadvantage of the acceleration potential as 
compared with the velocity-potential approach lies in the 
appearance of a more complicated kernel of the integral 
equation. 


INTEGRAL EQUATION 

For compressible flow, the integral equation as derived on 
the basis of the acceleration potential and given, for example, 
in reference 24 may be written as 


* ( 1 ) 

where w is the known vertical induced velocity or downwash, 
A y (£,i?) is the loading or pressure difference (positive Upward) 
at a point ({, 17 ) on the wing surface, V is the stream velocity, 
and K is the kernel of the integral equation. The pressure 
difference Ay (£, 77 ) is related to the bound vorticity through 
the equation 

Ap(Z,v)=pVT(Z, v ) 

so that equation ( 1 ) is sometimes written with re- 

2 

Op/> Yi) 

placed by — — The function K represents, physically. 


the contribution to the downwash at a field point (x,y) in the 
plane of the wing due to a unit pressure difference at a loading 
point (£, 17 ) also located in the plane of the wing. 

The kernel K is discussed in reference 24 and appears 
originally in an integral form as 


K\M, | (*-£), £( 2 ,- 77 )] 


—lim ^—5 e 

r-KI oz 2 


i~ [x-ivVx^+p’ j+ 7 ?u»] 

d 2 -«£(*-© n-i e vi> ' 


yx i +fl 2 (2/-7,) 2 +0V 


dX (2) 


where w is the circular frequency of the oscillating wing. 

In the form given by equation ( 2 ), the kernel function 
appears as an improper integral whose complete integration 
in closed form has not been accomplished. The kernel func- 
tion has been extensively treated in reference 24, however, 
and has been reduced to a form which can be evaluated. The 
reduced form of the kernel given in reference 24 contains 
Bessel and Struve functions and proper nonsingular integrals 
which can be handled by numerical means. The kernel con- 
tains singularities at y— 17=0 and z— £k 0 , but these singu- 
larities have been isolated and expressed in a form which can 
be handled in numerical procedures. Also given in reference 
24 is an expansion of the kernel into a series in powers of the 
reduced frequency, which for the present time serves as the 
most practical form for application of the kernel. 

The solution of equation ( 1 ) requires a determination of 
the unknown loading distribution Ay (£, 7 ?) subject to known 
boundary conditions. In view of the complicated nature of 
the kernel, an exact analytical solution of equation ( 1 ) does 
not seem possible. It appears necessary, therefore, to con- 
sider some numerical method of solution and several possible- 
approximate methods are discussed in the next section. 
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SOME APPROXIMATE PROCEDURES 

In searching for an approximate means of handling the 
integral equation, two approaches seem likely, one involving 
an approximation for the kernel and the other an approxima- 
tion for the loading. The first of these would involve a re- 
placement of the complicated kernel function by a simpler 
function which behaves in a similar manner (with regard to 
singularities, and so forth) but with which the integral equa- 
tion could be inverted and solved analytically for the un- 
known loading. Such a procedure has been used by Fettis 
(ref. 27) for the two-dimensional, subsonic, oscillatory case 
and, in part, by Lawrence (ref. 28) for the finite wings in in- 
compressible steady flow and by Lawrence and Gerber (ref. 
29) for oscillating finite wings in incompressible flow. The 
application of such a procedure to the case of compressible 
flow would be quite difficult, but is perhaps worthy of further 
consideration. 

A second possible approach involves leaving the kernel 
unchanged and selecting an appropriate expression for the 
loading which leaves only certain unknown coefficients to be 
determined. By such a procedure the integral equation can 
be reduced to a sum of definite integrals and the unknown 
constants can be determined by collocation, that is, by forcing 
the solution to fit the known downwash conditions at a num- 
ber of points on the surface. Although difficulties arise in the 
integrations because of singularities in the kernel and the 
loading, they can be handled in numerical methods. Various 
forms of this procedure have been used quite successfully for 
oscillating two-dimensional wings and for finite wings in 
steady flow. 

A logical extension of this second approach, the choice of 
a replacement function for the loading, to the three-dimen- 
sional problem would involve finding a means of handling 
the definite integrals which arise. It would seem possible 
to divide the wing into many small areas and, with an 
assumed form of the loading on each area, integrate numeri- 
cally over each area to find the downwash at a point on the 
surface. Such a procedure would lead, however, to a system 
of equations in the unknown constants of at least the order 
of the number of areas so that lengthy computations would 
arise which would require computing equipment of great 
capacity. 

In searching for a less cumbersome method, it is natural 
to study the possibility of utilizing some of the concepts of 
procedures which have already been developed for the case 
of steady flow. Several methods for steady flow seem adapt- 
able to the unsteady problem. Among those to be considered 
are the procedures of Falkner (ref. 25), Multhopp (ref. 30), 
Weissinger (ref. 31), and Schlichting and Kahlert (ref. 26). 
A brief discussion of these methods in relation to the objec- 
tives of the present investigation is given subsequently. 

The methods developed for the steady case by Falkner 
and by Multhopp represent attempts at lifting-surface pro- 
cedures for handling the integral equation. In both ap- 
proaches, the unknown loading is expressed in 'terms of a 
series of chosen modes of loading with unknown coefficients 
to be determined by collocation. In both methods, also, 
the double integrations are performed numerically and lead 


to sets of downwash factors. It appears that either method 
could be extended to the unsteady case and systematized 
through the development of suitable tables of the downwash 
factors. For the purpose of extending such procedures to 
the unsteady-flow problem, in the present investigation, the 
Falkner procedure has been considered preferable since con- 
siderably fewer tabulations are required. The development 
of a lifting-surface approach based on the method of Falkner 
is the main objective of this report. An extension of the 
steady-state procedure of Multhopp to the oscillatory case 
has been performed by Jordan and some results aro given in 
reference 32. 

It should be noted that W. P. Jones (refs. 2 and 18) has 
also treated the problem of determining the aerodynamic 
forces on wings of any plan form by using some of Falkner ’s 
concepts. The method differs from the one to be presented 
in this report and pursues a different numerical path since 
in the Jones procedure the basic integral equation was 
treated on the basis of the concept of the velocity potential. 
The use of the velocity potential results in a double integra- 
tion over the entire field of velocity discontinuity, that is, 
over the wing and its wake. The use of the acceleration 
potential, employed in this report, involves a double integra- 
tion only over the area of pressure discontinuity, that is, 
over the wing surface itself, since no pressure discontinuity 
exists in the wake. 

Other steady-state procedures mentioned previously wore 
those of Schlichting and Kahlert and of Weissinger. These 
methods represent lifting-line rather than lifting-surface 
approaches. In the Weissinger method a single lifting lino 
is placed at the X-chord position on the surface and the 
known downwash is satisfied at points on the %-chord 
position. The Schlichting procedure makes use of several 
lifting lines. The surface is divided into a number of span- 
wise strips of equal chord. A lifting line is placed at the 
%-chord position of each strip and the downwash conditions 
are satisfied at points on the /(-chord position of each strip. 
Since a lifting-line procedure might be expected to involve 
less labor than a lifting-surface approach, an attempt was 
made in the present study to extend the Schlichting method 
to the case of unsteady flow. The application of the Schlich- 
ting approach to the oscillating case is discussed in ap- 
pendix D. 

In the following two sections the lifting-surface mothod 
of the present report is developed, first for the case of 
subsonic flow and then for the limiting case of M=l. 

DESCRIPTION OF SURFACE-LOADING METHOD FOR 
SUBSONIC FLOW 

The solution of the integral equation for the case of an 
oscillating three-dimensional wing involves a double integra- 
tion over the plan form of the wing as indicated by equa- 
tion (1) : 

^^-=- JJ A V&v)K [AT, ^ (x-Q, £ (y~ v )j di- dr, 

S 

Of the three main ingredients of the problem — the downwash, 
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tko kernel, and the loading — which are related by this equa- 
tion, two can be considered as known. The kernel function K 
can be evaluated through the use of the forms given in refer- 
ence 24. The downwash w(x,y,t) is determined by the 
motion of the airfoil. This leaves as the only unknown in 
equation (1) the function Ap(£,ij), the magnitude of the load- 
ing at points (£, 17 ) of the surface. 

STATEMENT OF THE BOUNDARY CONDITIONS FOR THE DOWNWASH 

For linearized flow, the downwash function w(x,y,t) can be 
determined from the motion of the mean surface of the 
airfoil through the relation 

w(x,y,t) =~f(x,y,t) + V ^f(z,y,t) (3a) 

where, for harmonic motion, 

f{z,y,t)=J(x,y)e iut 

so that f(x,y) is the amplitude of the displacement of the 
mean surface from the equilibrium position. Once an appro- 
priate function for j(x,y) is chosen, the downwash is known 
and the loading can be determined from equation ( 1 ). For 
the usual modal type of flutter analysis, for example, f(x,y) 
could be found from the assumed vibration modes of the 
structure. 

Becauso the basic problem of determining the loading 
for a given downwash as expressed by equation ( 1 ) is linear, 
it is often convenient to consider the problem separated 
into a sum of individual problems, each associated with a 
particular selected type of motion or degree of freedom and, 
hence, with a particular downwash function. This implies 
that the total downwash can be written as a sum of individual 
downwash functions, so that equation (3a) appears as 

w(x,y,t)==(j^+V [Ji(x,y)+Mx,y)+ . . . +J q (x,y)]e iat 

=[wi(x,y)+w 2 (x,y)+ . . . +w q (x,y)]e fut (3b) 

where w q (x,y) is the downwash associated with the particular 
degree of freedom f q (x,y). The basic problem of this report 
' is then the determination of the loading Ap„(£,»j) associated 
with any one of the downwash factors w q {x,y). 

TREATMENT OF THE LOADING 

In addition to the downwash conditions set by equations 
(3), the loading Apfay) must satisfy various conditions at 
the edges of the plan form. In the chordwise direction, for 
subsonic flow, the loading must be infinite at the leading 
edge and vanish, in accordance with the Kutta condition, 
at the trailing edge. In the spanwise direction, the loading 
should become zero at the wing tips. 

Series expression for the loading. — As a first step in the 
determination of the unknown loading, it is assumed that the 
function A p(£,»j) can be expressed in terms of a series of both 
spanwise and chordwise pressure or loading modes, so chosen 
as to satisfy the edge conditions just discussed and contain- 
ing arbitrary coefficients to be determined. The spanwise 


modes are expressed in terms of the variable 17 . It is con- 
venient to express the chordwise pressure modes in terms of 
an often used variable 0 , related to £ and based on the 
circular cylinder enveloping each chord as diameter and 
shown by sketch 1 and in equation (4). 



£ —Jmi (v ) — 6 ( 17 ) cos 0 (4) 

where 

*(n)“! lMv)-Mv)] 

In terms of the variables 0 and 17 , the loading associated 
with a particular downwash distribution w q (x,y) may be 
expressed in terms of the summation of pressure modes and 
written as 

A p q (£, 17 ) Vs 2 — V* [cot | (a&+aS 1 77 +a&y i + . . .) + 

sin 0(af o +Gii’?+<4i’7 J + • • • )+sin 20((4>-|-(ij l i}-f 

• • • ) + • • -J (5a) 

where the superscript g identifies a particular downwash. 

The series form of the loading given by equation (5a) 
satisfies the Kutta condition (Ay=0 at the trailing edge) 
and has the desired type of singularity l[\x where x is 
measured from the leading edge. The terms in 17 are of such 
a form as to cause Ap(£,j 7 ) to become zero with infinite slope 
at the wing tips. The use of this series form for Ap(^,n) 
means the pressure distribution is essentially synthesized by a 
series of chosen pressure modes and the values of a nm deter- 
mine the contribution of each pressure mode to the final 
pressure distribution. 

Equation (5a) can be written in more concise form, as 
Ap<(M=!|K ± £ <4 n f n (Qg m ( v ) (5b) 
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... , /o(?)=COt| 

/*(£)= sin n8 
0m(v)=W«*— ? 

(For symmetrical motions of a wing, w a (x,y) =w,(x,—y) , 
the odd terms in 17 need not be retained. For unsym m etrical 
motions, w g (x,y) = —w ll (x,—y), the even-powered terms in -r) 
do not have to be retained. In the numerical examples to be 
given in a later section, symmetric motions are considered.) 

Illustrative expressions for the force and moment. — The 
specific problem is the determination of the coefficients a„ m 
in equation (5a) or (5b) and, as has been indicated, a par- 
ticular set of values for a„ m must be found for each down- 
wash w t (x,y). Once a set of values for a, m has been deter- 
mined, a particular loading A p g (%,y) is defined by equation 
(5a) or (5b). This expression for Ap„($,ij) gives the pressure 
distribution and can be integrated to obtain various desired 
section forces or total forces. For example, total lift and 
moment coefficients associated with the motion/, (x, y) can be 
found through the relations 

Cl. ^^JJa p, (I, y)d£ dy (7a) 

and s 

J(*-a)Ap,(*,„)#<fc, ' (7b) 

s 

where the moment is obtained about a spanwise axis £=a. 

By substituting equation (5a) into equation (7a) the 
following result for the lift coefficient is obtained which is 
valid for any plan form: 

—2 A 

^ A5= ~16~ ( 1 + 2 af, + 2 + «u) (8) 

The expression for the moment about the axis considered 
must be determined for each particular plan form. The 
moment coefficient will take on a form similar to equation 
(8) but with constants determined by the plan-form 
geometry. For a rectangular wing, for example, the moment 
coefficient for moment about the midchord (a=0) is 

u-M 

<r= — 22 " (16o«i+4flM+2a{i-|-8aa>+2oSj+a»i) (9) 

The equation for the moment coefficient for a delta-wing 
plan form is given in appendix C. 

REPLACEMENT OF THE INTEGRAL EQUATION BY A SUMMATION 
OF DEFINITE INTEGRALS 

In the rest of the analysis, it will be convenient to deal 
with a general downwash w(x,y) and to drop the subscript 
and superscript 7. In order to evaluate the coefficients a„ m 
for a particular downwash w„{x,y), the first step involves the 
substitution of the loading series (eq. (5a)) into the integral 


0=0) 
(n,p*0) f 


( 6 ) 


equation (eq. (1)) to give 

^% 2 ^ = 8^Xr' v/i¥=r ^[ COt l ( Ooo+Ooi’7+ao 2 u 2 + • • •) + 

8 

sin • • OOsin 20(o2o-f-fl2i??~f~ 

<w+ ...)+•• ^K^M.ytx-S^iy-ri^dl-dn 

( 10 ) 


or if the form for Ap({,ij) given by equation (5b) is used, the 
following is true: 


w(x,y) 

V 


s 

0,^-iy— dy 


(ID 


A significant step toward the desired solution of the in- 
tegral equation has now been arrived at because equation 
(10) or (11) is no longer an integral equation, that is, an 
unknown quantity no longer appears under an integral sign. 
Instead, the downwash w(x,y) is now given by a summation 
of definite integrals multiplied by various unknown constants 
a nm . These integrals are of the form 


JJ cot | V* V * 2 — pAf; ~ (r— £)-, ^r(y — >?)-ld£ dy 

s J 

JJsin nd y m -yJs 2 —y t K^M, ~ (x — £), y (y~ y) dy 


K (12) 


The numerical evaluation of these integrals is an impor- 
tant part of the labor of the present method, for once thou 
values have been calculated, one can proceed immediately 
to the determination of a nm and, therefore, to the determi- 
nation of the pressure distribution. In the method under 
consideration, the required integrations are performed nu- 
merically. Because the chordwiso and spanwise integrations 
are treated differently, they are discussed separately in the 
next two sections, the chordwiso integration being discussed 
first. 


TREATMENT OF THE INTEGRATION IN THE CHORDWISE DIRECTION 

For the purpose of performing the chordwise integration 
it is convenient to replace the right side of the integral equa- 
tion by a sum of separate integrals, each involving a par- 
ticular chordwise pressure mode and in which the spanwiso 
(77) and chordwise (£) integrations are separated. For this 
purpose equation (11) may be ■written as 

(13) 

where 

7o=J^/o©^[m, ^ (x~Q,f r (y-„)] d* (14a) 

and 

4=J^/.(f)^[M, y (s-{),f G/-1?)] dt (14b) 
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Remarks on methods of performing the chordwise inte- 
grations, — One could at this point evaluate the integrands 
of I 0 and I„ at small intervals of £ and integrate by numerical 
means. This procedure would be difficult and tedious, how- 
ever, since the kernel contains singularities and involves a 
large number of parameters and since the integrations would 
have to be earned out for many values of y. It is desirable 
to seek a simpler, but more approximate, approach which 
would be less laborious. 

In the present analysis, an approach developed by Falkner 
in reference 25 (and discussed more fully in ref. 33) for the 
steady case is used which greatly reduces the amount of 
laboj\ In essence it is assumed that in the integrals J 0 

Q 

and the loading modes /<>(£)= cot - and /„(£) =sin nd 

can ho replaced by an arbitrary but small number of re- 
placement loads which are chosen in order to satisfy some 
integrated properties of the continuous loading. The use 
of this assumption will result in the integrals J 0 and I n 
appearing finally as summations in the forms given by 
equations (20a) and (20b). 

The use of a set of individual loads in performing the 
chordwise integrations and the procedures for evaluating 
them introduce approximations which are somewhat arbi- 
trary. These approximations and their implications will 
be discussed later as they arise. They correspond to the 
assumptions made by Falkner and have been used in the 
present investigation for practical reasons. With high- 
speed computing equipment available, some of the approxi- 
mations in the development which follows might be avoided. 
It can be stated, howevor, that, in the cases considered, the 
end results seem to bo satisfactory. 

Calculation of the replacement loads. — In calculating the 
values of the replacement loads, it is first required that the 
sum of the individual loads associated with a particular 
chordwise pressure mode must equal the integrated value of 
the continuous chordwise loading. If the continuous dis- 
tribution is replaced by j replacement loads and 0 , o and 
O rn denote loads associated, respectively, with the modes 
cot 0/2 and sin 720, this equivalent load condition leads to 
the following set of equations: 

~^2G^^= f cot ^ sin 0 d0=7r (15a) 

>■■=1 Jo * 

and 

G rn = f sin nd sin 0 d£=~ (t 2 -= 1 ) 

r-l Jo * 

^30, n = f sin 710 sin 0 d£=0 (t 2 .)> 1 ) 

1 Jo 

The use of equations (15) imposes one condition on the 

replacement loads G. However, equations (15) contain j 
unknowns; therefore, at least j — 1 additional equations in 
terms of G are needed. In the present method the addi- 
tional equations are obtained by requiring the downwash 
produced at selected points on the chord by the individual 
loads to equal the downwash due to the continuous load 
distribution. In imposing this condition the two-dimen- 
sional kernel function K is used instead of the more compli- 


cated three-dimensional downwash factor defined by equation 
(2). The relation between the two- and three-dimensional 
kernel functions is given by the following equation: 

k[m, £ (*-£)]= “ (x-t), £ (7/-,)] dy 

where R [m,£(z-£)] is the two-dimensional kernel. 

The use of the two-dimensional downwash factor at. this 
point is an intermediate step used only for evaluating 0. 
In satisfying the actual downwash on the wing in the next 
section the calculated replacement loads are used •with the 
three-dimensional kernel function. This intermediate use 
of the two-dimensional downwash factor implies that the 
areas close to a point contribute the major portion of the 
downwash at the point. 

The,)— T additional equations in terms of G can be written 
by selecting j — 1 positions, or control points, on the wing 
chord denoted by x g at which to equate the downwash 
produced by the individual loads to that produced by the 
continuous loading. Although the lift conditions (eqs. 

(15) ) did not require any assumption regarding the location 
of the replacement loads, for the purpose of employing the 
downwash .condition such an assumption is necessary. In 
selecting this location an interpretation of the meaning of 
the replacement" -loads- -0 is somewhat-arbitrary: — The 
replacement loads O can be considered as loads distributed 
over a certain portion of the -wing chord, or they may bo 
considered as concentrated loads operating at particular 
points on the chord. If the concept of a load distributed 
over an area is taken, then the distance to the downwash 
point is physically identifiable as some average distance, 
which for this particular case is. moasured from the center of 
the load area. If the concept of the concentrated load is 
used, then the distance is measured from the point of appli- 
cation of the load. In either case, a reference point x, is 
selected as the point of application of the load G so that the 
distance to the downwash point becomes x g — x,. This 
distance is then used in developing the j — 1 downwash 
equations for the determination of G. 

For the cot 0/2 term of the loading series, the downwash 
equations are 

) GraK„— f cot sin 0 Kg d8 (16a) 

>■-1 Jo * 

and, for the loading sin t 20 , the downwash equations become 

QmKp— f sin 720 sin 0 dd (16b) 

•'—l J 0 

where one such equation is obtained for each control-point 
location x g selected. Tbe functions K t , and K& in equations 

(16) are defined as follows: 

Kc,=K(M,Z,) 

and 

K l( =K(M,Zi) 
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and represent the two-dimensional downwash factor tabu- 
lated by Schwarz (ref. 34). In these functions Z,=k(x t —x,) 

and Z[=k(x s —£) where k—j? and where a semichord b has 


been chosen for convenience as a reference length. 

It should be noted that the integrands in equations (16a) 
and (16b) become infinite at £=x and the integrals must 
be carefully treated. The method of handling these integrals 
is given in appendix B. 

Other conditions which would not bring in consideration 
of the two-dimensional kernel could presumably be used- for 
determining values of Q. For example, in place of using the 
downwash conditions, j — 1 additional equations could be 
written for each chord wise term of the loading series in which 
the first and higher (through j—l) moments of each load 6 
could be equated to the corresponding moments produced 
by the continuous loading. By such a process, however, 
the loads 6 would not be a function of either frequency of 
oscillation or of Mach number. Through the use of the 
downwash conditions, some effects of both frequency and 
Mach number are included at this intermediate stage. 

Illustrative equations for the loading functions. — In apply- 
ing the method of this report, three chordwise pressure 
modes have been retained. The continuous distribution in 
each mode has been replaced by four individual replacement 
loads. 

In order to evaluate the replacement loads, the lead 
established in the steady case with regard to locations of 
control points and loads has been followed. The loads are 
assumed to act at the %-, %-, %-, and %-chord positions, and 
the downwash conditions are applied at control points 
midway between the load stations, that is at the and 

%-chord positions. Replacement loads based on these 
chordwise positions have been calculated for M— 0 , 0.5, and 
1.0 and are presented in table I. 

One set of replacement loads must be determined for each 
chordwise pressure mode. For the mode involving cot 0 / 2 , 
the values of G, n are determined from the following equations : 


C T 0 

G r io - l~ <?2o+ ^(rb £?«> = I cot— sinftd0=7r 
Jo « 

6 r io5"n4~ GiaK\i-\- G\ w ^i 3 + GioSh— j cot ^ sin 0 dd 

6 r ioR , 21+ $‘0^22+ ^30^23+ G U) Ku= r cot|sin 0 ^ 2 ?.d 0 

£?10^1 + 6WR32+ ^ 30 ^ 33 + GitJ&i t -f cot^sinSR^dsJ 


Y ( 17 ) 


where the functions K used here and in equations (18) and 
(19) are defined after equation (16b). 

The first of these equations imposes the lift conditions 
stated by equation (15a). The remaining equations apply 
the downwash condition of equation (16a) at each of three 
control points; the second equation, for example, states that 
at control point 1 (at x s =xi in eq. (16a)) the downwash 
produced by the four loads at av Dl , x, m2 , x, m3 , and x ,_ 4 must 
equal the downwash produced at control point 1 by the 
continuous cot 0/2 loading. 

A similar set of equations can be written for each chordwise 


pressure mode by using equations (15b) and (16b). For 
sin 0 the set of equations is 


6?ai+ f?si+ Gti —J sin 2 0 dd=P 

$ 11 -^ 11 + GiiKii-\- G n K n -(- GiJ£\i= ^ sin* 0 E\t dd 

G n En-\- G-nK a + Gi \En= ^ sin 2 0 E# dd 

G n E 3l + G il E 3i + G z 1^33+ GiiE 3i —{ r sin 2 oE 3i ddj 
and for sin 20, 

Gi2-\-Gn+Gti-\-Qii= ^ sin0sin20d0=O 


( 18 ) 


G\JE\\-\- GaEn-{- G 33 Eis~\~ QiJEu= ^ sin 0 sin 20 E\t 
G n K n + G li E n + GnKv+ Gi % m. u =£ sin 0 sin 20 .S# 
GnE^i-}- GnK Si -\- GizEss- GnE 3i — J' sin 0 sin 20 Kit dd 


dd 

dd 


7 ( 19 ) 


Use of the replacement loads G in the chordwise integra- 
tion. — Once the replacement loads G have been determined, 
by solving equations (17), (18), and (19), continuous integra- 
tions indicated by equations (14) can be replaced by summa- 
tions of the products of the loads G and the three-dimensional 
kernel function as follows: 


Io=&£ GioE^M, ^ (x-x,),f (y-y)j ( 20 a) 
and I n (eq. 14b)) becomes 

G'nK^M, £ (x-x,), £ ( 7 /— 77 )] ( 20 b) 


where the factor b arises through use of the substitution 
£= — £> cos 0 . In effect the continuous integrations have been 
represented by the replacement loads and a mean value of 
the kernel. These expressions may be substituted into 
equation (13) as follows: 


V r~l n-0 m-0 


£ (x-x,), £ (y—n)\ dr, ( 21 ) 


leaving only the spanwise integration to be performed. 

SPANWISE INTEGRATION 

For the purpose of discussing the spanwise integration, 
equation ( 21 ) may be rewritten in simpler form as follows: 


where 



dr, 

( 22 ) 


AP, ( 17 ) =4 anmG, n g m (y) (23a) 

n«=0 to— 0 
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or 


tion (25a) can then be written as 


AP,(rj)=4pV 1 Tjs 3 —if [<? F o(aoo+aoii?+a 0 j’7 s )+<7»i(aio+ 

• • O - ! - • • •] (23b) 

Equation (23b) may be seen to correspond to the original 
form of the loading Ap(£,ij) given by equation (5a). In 
equations (23), however, the continuous chordwise-loading 
terms cot 0/2 and sin nd have been replaced by the loads 
6U and G, a . 

The integration of equation (22) may be handled by 
several procedures. A straightforward numerical integration 
could be performed for each value of x, by evaluating the 
kernel S' at a number of spanwisa stations. However, the 
kernel K contains some difficult singularities which have to 
be carefully considered. Moreover, in order to make tables 
of these integrals for general use, tables would have' to be 
made for every aspect ratio, sweep angle, Mach number, and 
reduced frequency. In order to facilitate the development 
of tables, it is desirable to make use of a mean-value integra- 
tion which involves integrating the kernel over a short 
segment of span and using the value of A P, at the midpoint 
of the segment. This method was followed in the present 
analysis and has the advantage that the integrals have to be 
tabulated only as functions of Mach number and reduced 
frequency. 

The downwash at a point (x,y) due to a small element of 
span of length 2e, the center of which is located at £=av and 
jj=u A r, is 

Aw(x,y) 1 p 5 

V ~8vg AP '(Vw F ’* i 24 ) 

where 

| (*-*,), ^ (y-„)] dy (25a) 

and where A P, (y N ) refers to the value of A P, at the mid- 
point of the span element over which the integration extends. 

The integral given by equation (25a) is of central impor- 
tance in the present method and represents a “downwash 
factor” giving the downwash at a point (x,y) due to a unit 
pressure loading acting over a span element of length 2e. 
The value of the integral depends only on the relative 
distance of the element from the point (x,y) and not on the 
spanwise location of the element on the wing. For con- 
venience of calculation, therefore, it is desirable to perform 
a coordinate transformation to the center of the element, 
so that the integral appears as a function of distances from 
the point (x,y) . 

In order to perform this transformation, let y'=r,—y N in 
equation (25a) where >7' is a new spanwise variable. Equa- 


^ rN= j_ c y (x—x,),^ (y— yt?— ■n') 1 dr,' 

=/_‘ ( K [M, ^ (*-*,),£ (2/'-V)] dr,' (25b) 

where y / =y—yN and where y' and r,' are shown in sketch 2. 



Sketch 2. 

SELECTION OF REFERENCE LENGTHS 


It is convenient at this point to choose certain reference 
lengths in both the downwash factor F rN , defined by equation 
(25b), and in the loading A P r {y N ), defined by equations (23). 
In the case of the loading, the wing semispan s is chosen as 
a convenient reference length. If a new variable n=y N ls is 
introduced, the loading A P,(y N ) can be written as 

AP,(/i) =4:8 pV* X) X). arnng m (jd) Q,n (26) 

n=0 ct«=*0 

where 

g* Ot)=M”Vi— m* 

and where p is measured to the midpoint of the spanwise 
segment. ■ _ 

In the case of the downwash factor F tN , the length e is 
taken as a reference length. The variables x, x,, y' , and r,' 
are then considered in a new sense to mean nondimensional 
quantities obtained by dividing the dimensional distances 
by the length e. The downwash factor F rtT then appears os 

K [MXx-x,) ,k(y' -V)] dr,' 

where 


The use of « as a reference distance introduces a factor 1/e 2 in the kernel (eq. (2)) as follows: 


K[MMx-£)My-v)]=ji- 


rx-t £ h-M-A'-WP (J,-,)HA>» 5 1 

lira rr — 5 I e ....... ... d\ 

ZH> 0 J— VX 2 +/3%-’?) s +/3 2 z 2 


(27) 


Equation (27) can be used to define a nondimensional kernel function K by writing 

K[MMx-x,)My'-v')]=kK[M; l(x-x,),: Hy'-n')] 


( 28 ) 
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so that K in equation (28) corresponds to the quantity in 
braces in equation (27) . ' _ 

With the use of equation (28), the downwash factor F, N 
appears as 

K [M,k(x—x r ) My'—n)] At 

Finally, a nondimensional downwash factor F tN can be de- 
fined as 

£ K [M, k(x-x r )My'-v')} An' (29) 

FINAL EXPRESSION FOR THE INTEGRAL EQUATION 

Equations (26) and (29) can now be substituted into 
equation (24) to write a final expression for calculating the 
downwash at a point (x,y) on the wing. If the downwash 
due to all the span elements over the wing is summed, the 
total downwash becomes 

£ S AP,(ji)F rif (30a) 

V 1 | fjf 

By using the expanded form for AP,(p) given by equation 
(26), the total downwash can be written as (for «=s/20) 

SSSE a m G,ngMF,„ (30b) 

V ATf n m 9 Vjt 


The only unknowns in this equation are the coefficients 
a M . By selecting a number of control-point locations (x,y) 
at least equal to the number of o^'s desired, equation (30b) 
can be written for each control point. Thus, a set of simul- 
taneous equations is obtained in terms of the coefficients 
a-am which can then be calculated. 

The next section is concerned with the calculation of tho 
downwash factors F rN - which have been defined in equation 
(29). 

DETERMINATION OF DOWNWASH FACTOR F . 

The value of the downwash factor F (eq. (29)) for the oscil- 
latory, compressible case is not available in the literature 
and consequently must be derived. However, the kernel K 
has been discussed in detail in reference 24, in which several 
forms of K are given. One form, as given by equation (20) 
of reference 24, would require a numerical integration in 
order to obtain F, since K cannot be integrated in closed 
form. Another form of the kernel, an expansion in terms 
of the frequency parameter k, is given by equation (64) of 
reference 24 and is the form employed in the present analysis. 
In reference 24, a discussion is given of the accuracy of the 
frequency expansion. It is indicated that for moderate 
values of k and*M the series is quite good for small distances 
x—Xy and y'—y'. As k and the distances increase, the 
series deteriorates. Overall accuracy is not seriously af- 
fected, however, since the major contribution to the down- 
wash at a point is from the nearby points on the surface. 


Accordingly, the series is integrated term by term to obtain an approximate value for F. Several of tho integrals aro, 
however, improper, due to the existence of singularities in the integrand, and the concopt of the principal value or finite part 
must be utilized. (See, for example, ref. 35 for a discussion of the finite part of integrals.) The results of the integration 
for an expansion to the 5th power of k is given in equation (31). (In presenting the expression for the downwash_factor F, 
it is convenient to drop the primes on the quantities x' and y' and the bar on the reduced-frequency parameter k. In the 
following expression, therefore, x, y, and k will be considered to denote the dimensionless quantities referred to « unit length.) 


FVv=J K(%o,y—y)di) 

3MW - log [(y-l)2? 2 -(y+l)I?,] y -32MV-fiW(5M 3 — 3) (-V+1)- 


,+ 


20 


’ 0(yf i)-M?i 
Xo+F 


906 ^. i-\3 i 0o . Zp-f-tfa r 0 +Ai 2M 3 

_0(y 1) l0 " (M+l)(y-l) +ZtS{y+L) i0g (M+iXy+i) 0 1 + * 1Og 0(y+l)+A 


2 M 1 


0(y — l)-(-i?2 


(2 — 4-M* — 4Af 4 ) x„ [(1/— l)f? 2 — (y-t-DtfJ+0 8 [(y+D 3 -(y-l) 3 ] (-27+|+|-7ri)+20«(j/-l)Mog^~^- 
20»(y+l) 3 log £*±2} +J^{ (3M’+4K< log g=§J§-8M«0*( 3tf+l)x 0 - 

^ [(y-l) 3 I? a -(y+l) 3 i? 1 ]+| (6 — 20 M 2 — 15 Af 4 ) [(y— 1)1?.— (jf+D^ffJ x 0 2 +^ (ZM 3 -1) [(y— 1)W— 

(y + 1) I?! 3 ] - 160*2/ (v*+ 1) +20 8 [(»+ 1) 1 (y- «<] } ) 


( 31 ) 
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where 

x 0 =x-Xr 7=0.5772157 

where x 0 =x—z, is the distance from the line of integration 
to a control point, made nondimensional by dividing by e, 
and is positive for locations of the control point behind the 
line of integration. The distance y is also nondimensionalized 
by e and is positive to the right apd is measured from the 
center of the integration segment as previously discussed. 
The reduced frequency k must also be based on the span seg- 
ment £ and is accordingly 



Of course the number of spanwise segments N can be 
arbitrarily chosen. In general, the more segments taken the 
more accurate the result. Falkner (ref. 25) has used e=s/ 20 
for most cases. As an example, the layout for a rectangular 
wing as used in the numerical example in appendix A is shown 
in figure 1 . 

The next section of the report is concerned with the appli- 
cation of the surface-loading method to the special case 
ilf=1.0. 

THE SURFACE-LOADING METHOD FOR SONIC FLOW 

The problem of calculating the forces for the limiting case 
of a wing which is oscillating in sonic flow is in general the 
same as for a wing at subsonic speeds with two exceptions. 
First, for the case of a sonic trailing edge, that is, when the 
trailing edge is perpendicular to the flow, the satisfaction of 
the Kutta condition (Ap—0 at trailing edge) is no longer 
necessary and another form of the series expansion for the 
loading Ap is applicable. (For the subsonic-trailing-edge 
case, i.e., when the trailing edge is not perpendicular to the 
flow, the loading series would, presumably, be of the same 
form as already shown for the subsonic case as given in eqs. 
(5) .) The second difference between the sonic and subsonic 
cases lies in the form of the kernel K, which is used for the 
calculation of the downwash factor F as defined in equation 
(27). The modification of the kernel K to the limiting case 
M~ \ has been performed in reference 24, but the integration 



Figure 1. — Typical layout for lifting-surface method applied to a 
reotangular wing. 


of the kernel with respect to the spanwise variable t\ must still 
be performed and is presented in another part of this section. 

The present method is based on the usual assumptions of 
linearized theory and the use of this approximate theory may 
be open to question. However, for the very thin wings now 
being used on aircraft, it is felt that the first-order effects as 
given by linear theory constitute the major effects for unsep- 
arated flow and adequate solutions will be obtained. 4 

For the case of a sonic trailing edge the form selected for 


the loading is 




q b^ 8 1,2 \jfj~ • • 

•) + 



• • ) + (£— fle) ii (Cho+ 


Chiy^r<te2i?-\- • • O'h • • • 



(32) 

or- 

^=||~SS 

g 0 Ln-o 0 

where 

/j«)2ra(’?)^j 


(33) 

/»(£-/*.)=(£-/,.)"■* 

and 

(n== 0, 1, 2, . 

• •) 





This series contains some aspects of the supersonic case 
as well as the subsonic. For instance, as in the supersonic 
ease, the loading is not zero at the trailing edge; whereas, 
as in the subsonic case, the loading becomes infinite at the 
leading edge. 

TREATMENT OF INTEGRATION IN CHORD WISE DIRECTION 

By following the scheme used for M<Cl, a set of simul- 
taneous equations for the chordwise replacement loads is 
obtained corresponding to equations (15) and (16). 

In calculating the values of the replacement loads, it is 
convenient to choose tbe semichord at some spanwise station 
as a reference length, as was done for il/<l, and to perform 
a coordinate transformation such that the origin lies on the ' 
wing midline. By such a transformation a new chordwise 
loading function /„(£) can be defined as 

/»(£)= (H-S)"^ (n=0,l,2, . . . ) ' 

where £ is now considered as a dimensionless variable based 
on a reference semichord and measured from the wing 
midline. In terms of this new variable f, the equations 
expressing the loading conditions (corresponding to eqs. (15) 
in the subsonic case) appear as 

£&»= f (l+fl"^# (34) 

>-■1 J —1 

5 Some conditions which most be met In tbe linearization of the governing equations in the 
neighborhood of Af-I have been presented by Miles (ret 38) and MoHo-Christensen (ref. 37). 

As stated by Miles, it Is required that the conditions 5, MS, KS, KMS«1 be satisfied and, 
in addition, that one or more of the following conditions be met: 

Uv-i|»5vj A»ji/j 

where 5 Is the thickness ratio, M is the Mach number, k is the reduced frequency, and A Is 
the aspect ratio. 
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and for the downwash condition (corresponding to eqs. (16) 
in the subsonic case) 


p'd+r-^ 

J — 1 


(35) 


A set of simultaneous equations as given by equation (36) 
must be formed for each value of n and solved for the associ- 
ated values of 6. 

As in the case of subsonic flow, the integrands of the inte- 
grals of equation (36) become infinite at £=x t . The evalua- 
tion of these integrals is discussed in the latter part of ap- 
pendix B. 

At sonic conditions, a disturbance cannot be propagated 
upstream, so that the labor of computing the loading func- 
tions is reduced as compared with the case for M<^ 1 since 
the downwash factors do not have to be determined and 
summed for conditions where the control point is upstream 
of the integration area. 

As an example, if the loading functions are assumed to 
act at the %-, and-%-chord positions and the control 

points are located at the X-, K-, and %-chord positions, then 

Kit=Sii=Kn=K 33 =Kn=Ki\=0 

since the terms represent downwash factors for control 
points ahead of the loading stations. Thus, the set of 
equations becomes triangular and can be solved by successive 
substitutions. 

TREATMENT IN SPAN WISE DIRECTION 

The integration in the span direction is carried out in the 
same manner as was done for M< 1, that is, the wing is divid- 
ed into many small segments, and, with the load A p assumed 
to be constant across the segment, the integration is perform- 
ed. The kernel for this case is given by equations (47) in 
reference 24. However, it is not possible to integrate this 
expression for the kernel in closed form and, as was done for the case when M< O, the kernel was expanded in powers of 
the frequency parameter k and integrated term by term. The concept of the principal part of a finite integral was used. 
The form of the integrated series is given in the following equation: 

Fu-i^j K[M,kx 0 ,k{y'—v')\M-idv' 


where the functions K v and Kg designate 
Kgr=K [M,k{x t —Xr)] a,_i 

and 

and are the two-dimensional kernel as defined by equation 
(B23) of reference 24 and given in appendix B by equation 
(B9). Unlike the two-dimensional kernel for the subsonic 
case, these functions are not tabulated for M— 1 and hence 
must be calculated. 

Tn expanded form and for the nth chordwise-loading term, 
the set of simultaneous equations is 

<?lB+6 r 2B+6 r 3n+<?I*=J' /n(£) 

6 f ln^21+6 r 2a J S21+G t 3B^23+6 , 4»^24— J d£ 


y (36) 


(y+ 1) log (rfl)M-( 2 /-l) log ( 2 /-l) s + 2 (y+ l)- 2 (y— l)]+^{3r [(y+l) 3 -(y-l) 3 ]} + 

(3 2y iti) [(y+l) 3 7 (y— 1 )^— 6r/+^(^-^i) | ^+ 1 ) t -jy- 1 ) S + 2[( y +l)3-( y --l)»3 log^- 

2[(y+l) 3 log (y+l) 2 — (y— l) 3 log (y-l) s -| (y+l) 3 +| (y-l) 3 l}+^{-§|+-i| [(y+l) 3 -^-!) 3 ]* 

(y+D 6 -(y-D 6 (y+i) 7 — (y— l) 7 , a. 8 ( r (37) 

480x o 8,0642V 3 ^1,920 \y-l y+l) J J 


where the parameters x 0 , y, k, and y are identical to those employed in equation (31). 
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REMARKS PERTINENT TO THE SURFACE-LOADING METHOD 

In previous sections of this report, some of the approxi- 
mations involved in the present lifting-surface method have 
.been discussed. It has been pointed out that certain arbi- 
trary features arise in perfor min g the numerical evaluations 
of tho integrals, particularly with regard to the chordwise 
integrations. The assumptions that have been made in this 
connection include the replacement of the continuous chord- 
wise loading by a number of individual replacement loads, 
the use of the two-dimensional kernel function in evaluating 
these loads, and the choice of a set of load stations at which 
the replacement loads are assumed to act. In addition, a 
set of control points must be selected at which the governing 
downwash conditions are to be satisfied. In the actual 
application of tho procedure these features give rise to certain 
problems which are discussed in the following.parographs. 

CONTROL-POINT AND LOAD-STATION LOCATION 

No attempt has been made in this report to determine the 
most favorable location of the control points or the optimum 
location for the stations at which the replacement loads are 
assumed to act. For such a study, systematic tables of the 
downwash factor are desirable and when these tables become 
available, it will be much easier to evaluate this aspect of the 
problem. In the present study, such tables were not avail- 
able, and the positions of the control points and the means of 
distributing the loading selected were the same as have been 
used for the steady case. It is felt that, at least for the simple 
rigid modes considered in the present study, the location of 
the control points is not critical. For more complicated 
modes of deformation other locations and additional control 
points might bo necessary. This question is worthy of 
further study. 

EFFECT OF TAPER 

In the calculation of the distribution of the loading across 
the chord, it is necessary to select a value of the reduced 
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k 

(a) Magnitude. 

Figure 2. — Variation of lift and moment with reduced frequency 

A=-2' 
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frequency to be used in calculating the two-dimen- 

sional downwash factors Kg, and Kgt in the integrals in 
equations (16). For a tapered wing, the problem arises as 
to what span position to use to obtain a reference chord for 
calculating k. 

In relation to this effect of taper, as a part of the present 
investigation, two separate calculations were made for the 
case of a 45° delta wing. In one case, the root chord was 
used as the reference chord, whereas, in the other, the chord 
at the midsemispan was used. It was found that the final 
results obtained for the pressure distribution and, conse- 
quently, for the lift and moment from the two calculations 
were almost identical. On the basis of this one test case, it 
may be inferred that the location of the reference chord is 
not an important factor. In the remainder of the delta- 
wing calculations presented in this report, the chord at mid- 
semispan was used. 

DISCUSSION OF SOME APPLICATIONS OF THE METHOD 

This section is concerned with a discussion of the results 
of calculations which were based on the lifting-surface method 
discussed. Calculations have been made for both rectangu- 
lar and triangular wings and comparisons of the results 
have been made with existing theoretical and experimental 
results where possible. 

rectangular wing 

In order to furnish a basis for comparison of results of the 
method with existing theory, calculations have been per- 
formed at M= 0 for a rectangular wing with an aspect ratio 
of 2 pitching about the midchord for various values of the 
reduced frequency k. The results are compared with the 
results given by Lawrence and Gerber (ref. 29) and are 
shown plotted in figure 2. In figure 2 (a) the magnitudes 
are plotted, and in figure 2 (b) the corresponding phase angles 
are plotted. Excellent agreement is obtained for both the 



k 

(b) Phase angle. 

fc] for a rectangular wing osoillating in pitoh about its midohord. 
M= 0. 
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magnitude and the phase angle. Although both the present 
method and the method of Lawrence and Gerber are approx- 
imate, the good agreement between the results promotes a 
feeling of confidence in both methods. It should be noted 
that the two methods are not similar and contain entirely 
different approximations. 

To obtain some effects of Mach number, calculations were 
made for the same rectangular wing with aspect ratio of 2 
pitching about the midchord for a range of Mach numbers 
at a constant value of k= 0.22. Results are shown in figure 
3. Included in the figure are the results of two-dimensional 
calculations. The magnitudes of the lift and moment are 
given in figure 3 (a) and the corresponding phase angles are 
given in figure 3 (b) . Calculations for three-dimensional flow 
up to and including M— 1 were made by the use of the 
surface-loading method. The results at supersonic speeds 
were obtained from reference 12. The variation with Mach 
number for both the lift and moment is approximately that 
which would be predicted by use of the factor ^/l —M 1 up 
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(b) Phase angle. 

Figure 3. — Variation of lift and moment with Mach number for a 
rectangular wing oscillating in pitoh about its midchord. A=2; 
fc=0.22. 


to M= 0.7. Both the moment and lift increase at M— 1 
and then drop off again at supersonic speed. Note that the 
two-dimensional lift coefficient has the same shape for the 
range up to M= 0.7 and that application of the aspect-ratio 


A. 1 

correction factor would apply fairly well. 


No such 


simple factor exists for the phase angles, and it would not 
be possible to correct two-dimensional results for finite span 
effects. 


TRIANGULAR WING 


Figures 4 and 5 present results for a delta wing with 
an aspect ratio of 4 oscillating in pitch about the midchord. 
In figure 4 is shown the lift and associated phase angles 
plotted against the reduced frequency k for M= 0. Results 
of the present analysis as shown by the solid line and those 
of Lawrence and Gerber (ref. 29) as shown by the dashed 
line are compared with some experimental results (indicated 
by the circles) obtained by Sumner A. Leadbetter and 
Sherman A. Clevenson at the Langley Aeronautical Labor- 
atory. It is noteworthy that the results of the two methods 
agree rather well, even with respect to phase angle. 

In figure 5 a corresponding comparison for the moment 
and its associated phase angle (for the same wing and for 
the same conditions as in fig. 4) is shown plotted against 
the reduced frequency k. Results of both analyses are in 
substantial agreement with respect to the magnitude of the 
moment. With regard to the phase angles, a significant 
difference between the results of the two theories occurs, 
although both theories indicate the same trend. The source 
of the difference can not be explained at the present time and 
will have to be resolved by further calculations and experi- 
ments. 
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wingjosoillating in pitoh about its midohord. A= 4; M= 0. 
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RESULTS OP MULTIPLE-LINE METHOD 

Figure 6 shows results of calculations for a rectangular 
wing with an aspect ratio of 2 at M= 0 based on the multiple- 
lino method described in appendix D. Results of the line 
method approach those of the surface-loading method when 
a fairly large number of lifting lines and control points are 
used. These results are discussed more fully in appendix D. 

CONCLUDING REMARKS 

The purpose of this report has been to present and describe 



k 

Figure 6. — Variation of moment with reduced frequency k for a 
delta wing oaolllating in pitoh about its midohord. A= 4; M=0. 


in detail a method for calculating the loading on a wing of any 
practical plan form which is oscillating in a subsonic or sonic 
stream. The method is presented in general form and some 
results of application are discussed. A sample calculation 
for a specific case is presented in an appendix. The method 
may be used for calculating the loading on elastic wings as 
well as on rigid wings. This feature makes the method 
adaptable to flutter calculations since it is possible to calcu- 
late the loading for the various modes usually assumed for a 
normal type of modal flutter analysis. The method can also 
be applied in principle, at least, to the combined aero- 
dynamic and structural problem in which the flutter char- 
acteristics are obtained directly by the use of structural and 
aerodynamic influence coefficients. 

The procedure used is based on linearized theory in which 
the usual assumptions of linearized flow, such as small 
thickness ratio of the wing, inviscid fluid, and so forth, are 
necessary. 

The method has been found to give good agreement with 
existing theory and experiment for low subsonic Mach num- 
bers for both rectangular and triangular wings. Results for 
high subsonic Mach numbers and a Mach number of 1 fit in 
well with theoretical results for supersonic flow but require 
further verification by comparison with experiment. 

In addition to the presentation of the lifting-surface 
method, a multiple-line approach is included in an appendix. 
Results of the line method approach those of the surface- 
loading method when a fairly large number of lifting lines 
and control points are used. 

It is realized that several variants of the procedure may be 
made and may be desirable for routine or systematic calcula- 
tions, particularly in view of the constantly increasing 
capabilities of automatic computing equipment. 

Langley Aebonautical Laboratory, 

National Advisoby Committee fob Aebonautics, 
Langley Field, Va., March 14, 1956. 



(a) Magnitude. (b) Phase angle. 

Figure 6. — Variation of lift and moment with reduced frequenoy k for a rectangular wing oscillating in pitch about its midchord. 
Comparison of results of multiple-line method with results of surface-loading method. A= 2; M =0. 






APPENDIX A 


NUMERICAL EXAMPLE OF THE LIFTING-SURFACE METHOD FOR A RECTANGULAR WING 

(A= 2; k= 0.22; M= 0.5) 


As an illustration of the calculation procedure of the lifting-surface method, the details of the procedure are presented 
in this appendix for a rectangular wing with A= 2 and £= 0.22, which is oscillating as a rigid wing about the midchord line 
in a stream flow of Af=0.5. For the case of subsonic flow the main equation to be dealt with is equation (30a) or (30b). 
The steps in performing the calculations are as follows: 

(1) The first step in the process is to divide the wing into a number of areas as shown in figure 1, to select the number 
and location of the control points, and, consequently, to determine the number of terms of the loading series (oqs. (5)) which 
are to be retained. For the present case, the wing was divided into four equal chordwise areas for the calculation of the 
distribution of the load in the chord direction. For the spanwise integration, the wing was divided into 21 segments of 
which 19 are segments of span 2e and two (one at each tip) are of span e. The semispan e of a segment is thus equal to one- 
twentieth the wing semispan. (In all of the calculations this quantity e is used as a nondimensionalizing factor so that the 
full wing span becomes 40 units. For a wing with A=2, as considered in this example, the chord is therefore 20 units and 
extends from — 10 at the leading edge to 10 at the trailing edge. The reduced-frequency parameter k (based on the half- 
chord) in terms of e becomes k=^=~ and for this example k= 0.22 and £=0.022.) 

Nine control points located at the and %-ehord positions and at three span stations 17=0.2, 0.5, and 0.8 of the 

wing semispan were selected. The selection of nine control points determines the minimum number of terms of the loading 
series. There can be more control points than terms of the loading series; however, a method of least squares would have 
to be employed for the final solution. 

In the present example, nine terms were retained in the loading series. Three chordwise terms, each modified by three 
spanwise terms, were used so that the series contains nine unknown coefficients and appears as 



(aoo+Ooi^+ffiM; 7 ? 4 ) +sin 0(aio+ai2i7 2 +anV) +sin 20{a 2o - J ra 2 A i -kchiV i ) 




(Al) 


It is only necessary to consider the even power terms in 17 since the loading is assumed to be symmetrical in the span direction 
about the midspan position. 

(2) The next step consists of calculating the replacement loads O from equations (17), (18), and (19). For the present 
case, there are three (n—0, 1, or 2) sets of four simultaneous equations as follows: 


(?!*+ f?2i»+ Q*n= SJniQ 


(—2.83500— 0.45000i) G ln - f (2.25500-0.47680i) 0 2a + (0.60680-0.29670i) G 3n + (0.29180-0.22086i) O in = fUZWii 
(- i,i7430-0.21760i) (? ln + (-2.83500-0.45000i) (? 2s + (2.25500-0.47680i) G 3n + (0.60680-0.296701) G in = ff n (£)K 2( d£ 
(-0.85000-0.09037i)G r lB +(-1.17430-0.21760i)6 f 2B +(-2.83500-0.45000^)G r 3,+ (2.25500-0.47680i)(7 4n =// B (f)^ )£ ^ 


r (A2) 


j 


The coefficients of the replacement loads G are values of the two-dimensional kernel K mn and have been obtainod 
directly from the table of reference 34. The kernel K mn is a function of M, and Z=k(x—Q. Consider as an example 

the center of area at % chord (%—■ and a control point at % chord (x=— 5). For £=0.022, Z= (0.022)^— 5+^^=0.066 

and for 111=0.5, the downwash factor H u = —2.83500— 0.45000i may be read from the table of reference 34. Similarly, 

for the same values of £ and M and for the center of area at % chord an( l the control point at midchord (a;=0), 

Z=0.022 ^0—^^=— 0.165 and, therefore, i£u=0.60680— 0.29670i. 
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The integrals on the right of equation (A2) are given as 
follows for the three sets of equations corresponding to 
n=0, 1, and 2. (See appendix B for discussion of method of 
integration.) 

0 

For 7i=0, /o(£) =eot - and 

nr q > 

J cot-sin 0d0=ir 

f "cot ^En sin 0 de= -2.88354—1.34125i 
Jo * 

- (AS) 

f cot sin 9dB=— 2.97504— 0.96137* 

Jo 2 

f" cot sin 0 d0=— 3.03074— 0.57320* 

Jo 2 

For 7!=l,/i(f)=sin 0 and 
J sin 2 6 d0=^ 

sin 2 0 E u d0=O. 589 12— 0.63606* 

I (A4) 

J^sin 2 0 Kn dQ— — 0.45095 — 0.67032* 

f "sin 2 0 K 3i dB=- 1.48547-0.55366* 

Jo J 

For 7*=2, / 2 (£)=sin 20 and 
J sin 20 sin 0 dd = 0 

J sin 20 sin 0 En d6= — 1.03882 — 0.20631* 

* (A5) 

J' sin 20 sin 0 E^ dd=— 2.02994+0.04145* 
f" sin 20 sin 0% f dd=— 1.02892+0.28861* 

J° ) 

The values of # as computed from the three sets of simul- 


taneous equations formed by substituting equations (A3), 
(A4), and (A5) into equation (A2) are 

6 10 = 1.72184+0.00224* 
#^,=0.73737-0.01130* 

(?3o=0.44052-0.00484i 
#40=0.24186+0.01389* 

#„=0.30752+0.02241* 
# 21 =0.47670+0.00582* 
#,,=0.47702—0.00583* 
#41=0.30956—0.02240* 

#13=0.45744+0.01103* 
#12=0.23804—0.01055* 

#33= -0.23753-0.01073* 

#44= -0.45795+0.01024* 

When the values of # are obtained, the loading function 
AP,(ji ) appearing in equation (30a) and defined by equation 
(26) may be written as 

AP,(p) =4:S~pV i -\/l — F [#ro(aoo+ao2*i 2 +n<MM 4 ) + 

#>i (am + «i id 1 + ttuM 4 ) + #+&2 o + 0 * 2 + + u ^h 1 ) ] (A7) 

where v=l, 2, 3, or 4 corresponding to the 4 areas in the 
chord direction. 

Exami nation of this equation shows that certain products 
of #raiU m Vl — 810 needed. These products correspond to 
the products g m (d)G n in equation (30b). It has been found 
convenient to arrange these products in a certain form for 
later calculations. This form is shown in table II where 
values of the products are given for the complete system of 
replacement loads. 

(3) The downwash factors F must now be determined. 
These factors may be calculated by the use of equation (31), 
which was used for this example. With the adaptation of 
high-speed computing machines to the problem, the kernel 
in the form of equation (20) of reference 24 may be numeri- 
cally integrated without, perhaps, the frequency and Mach 
number restriction of the series. However done, the distance 
between the centers of areas and the control points, x 0 and y', 
must first be determined. For the present case, the y' dis- 
tances for control point 1 and the first chordwise area corre- 
sponding to < 7 = 1 , for example, are 



The distance x=2.5 is constant since the wing is rectangular. 
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Since the loading in the example is symmetrical about the midspan, the loading factors which multiply F in equation 
(30b) are the same for equidistance on each side of the midspan. Therefore, the values of F which are located at equal 
distances from the midspan may be computed and added before multiplication by the loading; consequently, distances y' B 
and y' L have been given in the example where refers to the distance from a control point to the centers of the areas to the 
right of the wing midspan and y' L refers to distances to areas to the left of the midspan. The final form of the F factois is 
given in table HI where, for each entry, the two values' of F for y's and y' L have been added. 

(4) The next step is the multiplication of the loading factors of table II by the downwash factors F of table HI in accord- 
ance with equation (30b) which results in a 9 by 9 matrix of complex elements. This is given by the left side of the following 
equation and the downwash conditions to be discussed subsequently are shown on the right of the equation: 


(1.07713+0.181567)a,)o+(— 0.05279+0.095897)a,o+(0.33112+0.041297)a2o+(— 0.14365+0.043877)ao2+ 

(— 0.05511+0.01465f)a u +(— 0.02878+0.00716i)a 22 +(-0.07449+0.01803t)a w +(-0.02325+0.00588T)au+ 
(-0.01569+0.002621)054= -(0.025-0.002757)2x 

(1.20324+0.023397)000+ (0.27977+0.08613f)o 10 + (0.60908— 0.03442t)a5o+ (-0.21994+0.061731)005+ 
(—0.08525+0.022087)o 1 5+(— 0.03123+0.009217) O55+(-0.10617+0.027587)aoi+(-0.03939+0.009247)a 1 4+ 
(—0.01984+0.004347)054= — (0.025+07)2*- 


(1.28588— 0.123237)ooo+ (0.60852+ 0.031327) Oi 0 + (0.32497— 0.104817)aso+(— 0.26574 + 0.087927) Oo5+ 
(-0.11441+0.033047)o 1 j+(— 0.02668+0.011777)a55+(-0.12604+0.040137)a w +(-0.05509+0.014477)a 1 4+ 
(— 0.01473+0.006157)054= — (0.025+0.002757)2*- 


(0.90622+0. 110767) Ooo+ (0.01456+0. 066727) O! 0 + (0.27334+0. 0224 97) O5o+ (0.95139+0. 001487) a<) 2 + 
(0.11217 +0. 022727) O 15 + (0.25999 +0.003957) 055 + (0.65180 — 0.007157) Oo4+(0.08428 + 0.012287)a,4+ 
(0.17638+0. 001167) 054 = — (0.025 — 0.002757)2®- 

(1.11514 — 0. 019627)a w +(0.29257 + 0. 052947) a 10 + (0.47578 — 0.034107) 050+ (1.24921 — 0.130107) 0^5+ 
(0.37821 — 0.001327)045+ (0.41833 — 0. 042717) O53+ (0.86338 — 0. 096667) 004+ (0.26493 — 0. 004907) O14+ 
(0.28075-0.029777)054= - (0.025+07)25r 

(1.23341 — 0.156247)aoo+ (0.56690+0. 001267) Oio+ (0.26664 — 0.088067) 050 + (1.42203 — 0.281897)0^+ 
(0.63972 — 0.060637)Oia+(0.25104— 0.088387) 055 + (0.98606 — 0.201647)a w + (0.44244— 0.045957) O 14 + 
(0.17013 — 0.060287) Os 4 = — (0.025+0.002757)2x 


(AS) 


(1.01349 +0.157727)000+ (—0.03208+0.085737)040+ (0.30911+0.035207) O5o+ (0.25606+0.035427) OW+ 

(— 0.00455 +0.019797) Ois+ (0.07656+0.008157) Oss+(— 0.04884+0.023257) Oo4+(—0.02859 + 0.007977)oi4+ 
( — 0.00894+0. 004487)054= — (0.025 — 0.002757)2 tt 


(1. 15764+0. 009757) O00+ (0.27905+0. 074817) O10+ (0.56079 — 0.033817) O50+ (0.29204 — 0.001927)003+ 
(0.07234+0.016617)ou+(0.13788— 0.008797) 052+ (-0.09262+0.028617)oo4+(— 0.03866+0.011337)0,4+ 
(— 0.00479 +0.003937)024=— (0.025+07)2*- 


(1.24922 — 0.131267) Ooo+(0. 58649+0. 021877)o,o+ (0.30293 — 0.098267) Oto+ (0.31639 — 0. 037207) Oos+ 
(0.14828+0.003217)o^+(0.07495 — 0.024527)055+ (—0.11594 + 0.040077)00,+ (—0.04833+0.016037)0,4+ 
( — 0.00800 +0. 004027) Os«= - (0.025 +0.002757)2 tt 


(5) The boundary conditions w at the various control points are determined from the motion of the wing. For the 
present case, oscillations about the midchord line of a nondeforming wing are considered. 

Since f(x,y)=a(x— a) and, according to equation (3a), 


then (for downwash, lift, and displacement positive in the same sense) 

w ( x ,y) — — [F«+ a{x—a)\ 
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and for harmonic motion 


For the example considered the values of the a m are 


(A9) 

It is convenient to divide the simultaneous equations by the 
factor 40/2 v which appears in equation (30b). Therefore, 
the downwash conditions for rotation about midchord 
(a=0) are 


wfruVi) w(x 4 ,y^) _w(x 7 ,y 7 ) 


v« 

~ Fa " 

Fa 

w(Xi,y s ) 

w(x s ,y t ) 

w(x s ,y s ) 

Va “ 

” Vd 

Fa 

w(.V3,Vs) 

w(x 0 ,y e ) 

w(xv,y<>) 

Va 

~ Fa " 

Fa 


(0.025— 0.00275^)27^ 


( 0 . 025 + 0 ^) 2 •?^ 


r(AlO) 


(0.025+0.00275f)2ff 


Ooo= -0.150858+0.0334381' 
o 10 =0.045628— 0.129162i 
020 = — 0.000862 - 0.008337^ 
002 = -0.034965+0.007107-i 
0x2=0.057334— 0.031615i 
022 = 0.014221 - 0.004922i 
Oo4=— 0.025231+0.003860i 
Oh= 0.058109 — 0.006198i 
024=0.045049— 0.007729i . 


(All) 


(7) These coefficients may be used in equations (5) to 
determine the pressure distribution and, in turn, to obtain 
section forces, total forces, and so forth. The total lift and 
moment, for example, are given by equations (8) and (9) and 
for the sample case are 


where x ( and y t are the coordinates of the control points. 

(6) By the use of the boundary condition w/Vcc, the set of 
simultaneous equations (eqs. (A8)) may be solved for a m . 


10^1=2.632 


|Or.«|=l-594 


<h, a = 194.43° 
4>u,a= 349.14° 


(A12) 



APPENDIX B 

. TREATMENT OF CERTAIN INTEGRALS WHICH CONTAIN SINGULARITIES THAT ARISE IN THE CHORDWISE INTEGRATION 


In the treatment of the chordwise integration in the sur- 
face-loading method, certain integrals arise which contain 
singularities (eqs. (16)) and which must be given special 
treatment. In the determination of the coefficients 0, x 
appearing in these equations, it is necessary to evaluate 
integrals of the form 

In— f /*(£)-&+ d£ (Bl) 

U chord 

where /„(£) is the nth typical chordwise pressure mode in the 
series expression for the continuous pressure distribution and 
Ktf—K [M,k(x g — £)] is the two-dimensional kernel function. 
In this equation, the kernel function becomes infini te at 
£=r r so that the integrand is singular'. In the following 
sections, the methods of evaluating the integrals and taking 
care of the singularities are discussed. First the case of 
subsonic flow and then the case of sonic flow is discussed. 


CASE OF SUBSONIC FLOW 

For the case of subsonic flow, the integrals occurring in 
equations (16a) and (16b) are 

/ 0 = f cot | sin 0 dd (B2a) 

1,,=^ sin n6 sin 0 dB (n= 1,2,...) (B2b) 

Since the kernel 'in the integrands of / 0 and /„ becomes 
singular at £—x g) it is necessary to separate the kernel into 
singular and nonsingular parts. The singular part can be 
integrated in closed analytic form. The nonsingular part 
can be handled routinely and accurately by numerical meaos. 
The separation of the singularity has been accomplished by 
Schwarz (ref. 34) in the following manner: 

K(M,Z) F(M) +iQ{M) log, \Z\J+K l (M,Z) (B3) 
where 

Z=k{x t —£) 

The singularities now appear only in the quantity in brackets, 
and the nonsingular part K x has been tabulated by Schwarz 
(ref. 34). Substituting this expression into the integral of 
equation (B2a) gives (after sotting £=— cos 6 and x t =— cos 9 0 ) 

cot | sin dK(M,Z)dB= J (1+cos 0)K 1 (M,Z)dO+ 

£ (l+cos e ) [- log, |Z|] de (B4) 

The first integral on the right-hand side of equation (B4) 
must be evaluated numerically. The second integral may 
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be found analytically and has the following value : 

£ (1+cos 6) log, |Z|] dO=lF(M) + 

iO(M) [—■ a- log, 2—i r cos 0 ,+tt log, &] (B5) 

The treatment of the integral of equation (B2b) for n= 1 
proceeds in a similar manner to yield 

f sin 2 0K(M,Z)dd= f sin 2 0K x (M,Z)dB— r F(M) cos 0„+ 
Jo Jo * 

y Q(M)(— log, 2+i cos 20, + log, (B6) 

CoiTespondingly, for n=2, 

f sin 20 sin dK(M,Z)dd= f sin 20 sin 0 Ki(M,Z)dO— 

Jo Jo 

| F{M) cos 2 0 O + j G(M)(£^~co 5 0 o ) (B7) 


CASE OF SONIC FLOW 

For the case of sonic flow, the form of the series expression 
for the chordwise loading and also the form of the kernel 
function differ from those of the subsonic case. 

For M— 1, the set of equations to bo solved for the chord- 
wise loads (?,„ is given by equations (34) and (36). In 
equation (35), integrals of the following form appear: 


In=£‘a+k) n 2 K g & (n= 0 , 1 , 2 ...) 


(B8) 


AtM=l, the two-dimensional kernel function appearing in 
these equations is given by equation (B23) of reference 24 
and may be written as 


jr _2-/v 


(l+i)kg 2 

n/kx,, 


(1 


y 039 ) 


where 




and where and a] ' e Fresnel integrals 


defined by 


and 


C(a)=j cos ^ t*dt 
S(ot ) sin dt 
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As in the case of subsonic flow, the kernel function can be 
separated into singular and nonsingular parts. For this 
purpose, equation (B9) may be rewritten as 


where 


so that 




(BIO) 

(BID 


(1 - i )V^fe-“{c( A /^)+<s( A /S)]| (B12) 

The integral /„ (eq. (B8)) can then be rewritten as 


(B13) 


The first integral on the right of equation (B13) i 3 non- 
singular and can be evaluated numerically. The second inte- 
gral contains a singularity and can be integrated analytically. 

For 7 i= 0 , the second integral on the right of equation (Bl3) 
is found to be 


(■’■KeL dl = 0 r-- 

J — i VD+fl */— 




where 


.= Cir 


l V( 1+5) 

fVl a+i) * 

Correspondingly, for n= 1, the integral becomes 
For n=2, there is obtained 

cf x ‘ aa=c{x t +iy f 


£+t^-C{x t + l)f 


(B14) 

(B15) 

(B16) 

(B17) 


APPENDIX C 

CALCULATION OF MOMENT ON DELTA WING 


As mentioned in the text, the expression for the moment 
depends on the particular plan form. As an indication of 
the general method, the specific case of moment on a 45° 
delta wing is given in this appendix. The moment is calcu- 
lated about a line £=a as shown in sketch 3. 



£ 

Sketch 3. 

The expression for the moment as given by equation (7b) 
may be written as 

M a =lSb Cu, a = J'J (£— a)Ap a d£oh? (Cl) 

s 

When equation (4) is used (for a 45° delta wing) the dis- 
tance of an arbitrary point ( 77 , £) from the line £=a. in terms 
of the variable 6 is 

a — a=s [| tan A— 0r — y tan A^ cos 6— a -yj (C 2 ) 


Substituting Ap (from eq. (5a)) and equation (C2) into equa- 
tion (Cl) and noting that 

da=s (~~2 tan 0 dQ—b(v) sin 6 dd 

results in the following: 

Af a =8s 3 pF 2 J o -Jl—if [cot | (^+<W+<W) fl- 
am diaio+au^+auv*) fl- sin 2d{a K + (htrf+Ouv*) J[| tan A— 
(~~~ 2 tan A^ cos d—a sin d dd drj (C3) 

In this example, three ehordwise and three spanwise terms 
have been retained in the series for the loading Ap. 

Performing the indicated integration in equation (C3) re- 
sults in the following equation for the moment on a 45° delta 

wing [by making use of the fact that tan A=2 
M.= 8 rt. P FV [o»(i-|-«f)+a« Q-§- 

B n)+““ (sHan 0 §)+“» (r - ® 1 } 


Ul2 


(w»D 


fl-Ui 


+ 


4 (105 a 32)"^" 320 (12 ie) + 


022 (30 64) +ai4 (l05 12s)_ (C4 ' ) 


400104 - 08 - 


-04 
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APPENDIX D 

DESCRIPTION OF A MULTIPLE-LINE METHOD 


Since a lifting-line method of handling the integral equa- 
tion (eq. ( 1 )) may often he simpler than a lifting-surface 
procedure, it was considered interesting to investigate such 
a method and to compare results with those obtained by the 
surface approach. The purpose of this appendix is to pre- 
sent a lifting-line method based on a steady-state procedure 
developed by Schlichting and Kahlert (ref. 26) and to give 
some results obtained by it. 

DESCRIPTION OF THE METHOD 

In the Schlichting procedure for steady flow the lifting 
surface is divided into n spanwise strips of chord- c'=c/n, 
where c is the total chord. A lifting line is placed at the 
X-chord position (c'J 4) of each strip and the downwash is 
satisfied at control points on the %-chord position of each 
strip. This K-chord and X-chord location, respectively, of 
the lines and control points is an essential element of the 
procedure. It is a known fact that for a two-dimensional 
flat plate in steady incompressible flow the exact value of 
the lift can be obtained by placing a single lifting line at 
the K-chord position and satisfying the downwash at the 
X-chord position. Schlichting and Kahlert made use of this 
fact in developing a procedure for finite swept wings. 

In adapting Schlichting’s method to the oscillatory case 
the same placement of lifting lines and control points has 
been used and is shown in sketch 4, where, for example, the 
wing has been divided into two spanwise strips: 


assumed that the loading on the line l n can be represented 
by a series of the form 

A p M = pV 2 s(a„i sin 4>+a K3 sin 34>+a^ sin 6 $+ . . .) (Dl) 

which contains only one variable, the spanwise coordinate 4>. 
Since a separate loading function of this form is ■written for 
each lifting line, the double integration of equation ( 1 ) is 
reduced to a suin of single integrals. The integral equation 
then appears as 

^T^= 8^- 2 ? jspan *P*Kl M M*-Xr)MV-v)} dr, (D 2 ) 

where x, is the chordwise coordinate of the nth lifting fine, 
and where K, the non dimensional form of the kernel func- 
tion defined by equation (28), has been employed. 

When A p n in equation (D 2 ) is replaced by the series ex- 
pression in equation (Dl), the integral equation is reduced 
to a summation of definite integrals multiplied by the co- 
efficients. a nI , a * 3 , and so forth, as follows: 

S [a.i sin $ K(M, kx 0 , hj 0 ) dr, + 

®n3 /: sin 3<PK (M, kx 0 , hj 0 ) dv + • • -J (D3) 

where 

x 0 =x—x r y„=y— v 



In this example two lifting lines h and U axe considered 
that he, respectively, at X and % of the total chord. Two 
control points are shown, by aq at %c and at %c. In the 
sketch $ is an angular spanwise coordinate related to a nondi- 
mensional variable t) (referred to the semispan) byr?=— cos 4>. 

In the original form of the integral equation given by 
equation (1), the continuous pressure distribution Ap is re- 
placed by the sum of loadings on the individual lines. It is 


This equation sums the downwash at a particular control 
point due to all the lifting lines and the various types of 
loading, sin $, sin 3$, and so forth, on each line. The 
specific problem is the determination of the coefficients 
a„ u a * 3 , and so forth. In order to obtain the coefficients a, 
equation (D3) must be written for each of a number of con- 
trol points. This leads to a set of simultaneous equations 
which can be solved for the values of a. At each control 


point the function 


w(x,y) 

V 


is determined from the motion of 


the wing in the same manner as in the surface-loading method. 

Once the values of the coefficients a nX , a^, and so forth, 
have been found, they can be used in equation (Dl) to define 
the pressure distribution on any lifting line and in turn to 
give various force or moment coefficients. As an example, 
the total lift on a rectangular wing can be obtained from the 
relation 


C L 


n 

qS 


1 


~qS 



Ap, dy 


(D4) 


22 
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whore L„ is the lift on the nth lifting line. The total pitching 
moment about an axis x=a for a rectangular wing is given 

by 

Ma 
~gSb 


Cu=- 




0D5) 


where x„ is the chordwise coordinate of the nth lifting line. 

The main computational problem in the procedure is the 
evaluation of the spanwise integrals of equation (D3), and 
the handling of these integrals is discussed in the next 
section. 

EVALUATION OP THE SPANWISE INTEGRALS 

The spanwise integrations in equation (D3) are performed 
numerically and, since the kernel function has not been 
tabulated, it is necessary to make use of the series form of 
the kernel given by equations (31) and (54) of reference 24- 
In this form the nondimensional kernel function K appears 
as 


K(M,Jcx 0 ,ky 0 )=e <te » 



1 ^ 1 

\ yo\l^o 2 +fi% 3 J 



[ M +^ 2 (2 7 ) ^ log 2 ( 1 - 


irifi 2 


ixJ+Wyo 


-P 1 log (Vx 0 2 +/3V- 



0D6) 


where all distances are referred to the semispan s. 

For the integrals of equation (D3) which relate to lifting 
lines behind the control point (x o <0), the products of the 
loading modes and the kernel K as defined by equation (D4) 
can be evaluated at a number of values of y and the numerical 
integrations can be readily performed. When the lift i ng 
line lies ahead of the control point (r o >0), however, the 
kernel becomes infinite at y—y, and the singularities must 
be carefully treated. 

For the case r„>0 with y—y, singularities arise in the 
terms 


I Ntf+W +s<A 

V yoS^7+^v} ) 


5=—^ log (V*7+^V — : x °) 


(D7) 


Extraction of the singularities. — The singularities of the 
terms A and B can bo extracted and, when combined with 
the loading terms sin $, sin 3$, and so forth, are integrable 
in closed form. 

For this purpose, term A can be rewritten as 

a = r_/vgZjw+^\ + Ai_^ P8) 

L Wv^T+isv/ y°\ y° 


where the quantity in brackets is finite for all values of 
y and where 

' ‘ ^ 

2x 0 2 


lim 




The remaining term on the right of equation (D8) contains 
the singularity and will be handled analytically. 

The singularity of term B can also be isolated by writing 


B =—J log (Vav'+PV— x o)= J lo g( 


■yjxj+pyj+x, 


P, s 

-2 log y<? 


(DIO) 


where the bracketed quantity is finite for all values of y 
and where the remaining term can be handled analytically. 

Special form of the kernel for the case x o >0. — The ex- 
pressions for the terms A and B obtained in equations (D8) 
and (DlO) can be used to define a special form of the kernel 
function for use in the case r o ^>0. For this case the kernel 
can be written as 

K {M, kx 0 , bj 0 ) X ' >0 =K (Xo>0,y o ) +K' (?o>0,yo) (Dll) 
The term K(x„^>0,y o ) is nonsingular and is defined by 

K(.^o>0,yo) —K(M,Jcx 0 , Jcy 0 ) 

|[log +log ^° )] } (D12) 

where K(M,k,x 0 ,y 0 ) is the form of the kernel defined by 
equation (D8). 

The second term on the right-hand side of equation (Dll) 
contains the singularities and is defined by 

K'{x o >0,y o )=-e-**° Q+| log y 0 s ) (D13) 

The form of the kernel given by equation (Dll) is used in 
evaluating the integrals of equation (D3) for av>0. 

Performance of the spanwise integrations. — It is recalled 
that for the case x o <^0, the integrals of equation (D3) can be 
readily evaluated numerically by making use of the form of 
the kernel given by equation (D8). For the case r o )>0, the 
integrals of equation (D3) can be handled by using the special 
form of the kernel defined by equation (Dll). With this 
expression, a typical integral of equation (D3) for Xo)>0 can 
be written as 

J sin ra$ K(M,kx 0 ,ky<,)x.>Qdy= j sin m$K(x o ’>0,y o )dy+ 


/: 


sin m$^ / (r o >0,y o ) dy (D14) 
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where the functions K{x 0 ^>Q,y^ and IL'Ca^X),^) are defined, 
respectively, by equations (Dl 2 ) and (D13). The first 
integral on the right of equation (Dl4) can be readily evalu- 
ated numerically. The second integral contains the singu- 
larities and must be evaluated analytically. Its values 
have been determined for two pressure modes, sin <f>=V 1 — 17 1 2 
and sin 3$ = ( 417 s — 1 )VT~ V and are 

J* V 1 — t; 3 4 5 S' > 0 ,y 0 ) dn=— 2?r+ 

F-Vi —y 3 [(y+D log (y+i)— (y— l) log (y—i)—2]+ 

f (VW-VI=?) log Vo 2 dv } (D15) 

and 

J_VW ( 417 s — 1)2' {Xo'P’0,y o )d'n=—e- ttx °J^ — 6tt ^^° + 

^Vl —if (4i/— 1) [()/+ 1) log (y + 1) — (y — 1) log (y — 1) — 2] + 
j [Vw (4i, 2 -l)-^/I=y (4i/-l)] log (D 16) 

where 

4> 0 =cos -1 y 

For the special case of a control point at y— 0 and £<,>(), 
equation (D15) becomes 

J\/W K (Xo>0,-Ti)dy=-e-**° [-2ir-| (j+irlog 2 )] 

CD17) 

and equation (Dl6) becomes 

(4^-1) VW K'(x 0 >o,- v )dv=-e-^ [eir+l (J)] 

CD 18) 

For the case of y^O and r o >0, the integrals remaining on the 
right of equations (D15) and (Dl6) are finite and can be 
evaluated numerically. 
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TABLE I.— VALUES OF <?,„ 
[n=0, 1, 2; ,=1, 2, 3, 4] 

(a) Af=0 


k 

(?10 

0 M 

(?jo 

G, o 

0 

0.02 

1. 718194+0. 0006161 

0. 736332-0. 000638* 

0. 441765-0. 000412* 

0. 245302+0. 000434* 

0. 06 

1. 719093 + 0. 001479* 

0. 736393-0. 001990* - 

0. 441523 - 0. 001162* 

0. 244586+ 0. 001673* 

0. 10 

1. 720391 + 0. 001845i 

0. 736438-0.003446* 

0. 441159 - 0. 001816* 

0. 243605+0. 003417* 

0.20 

1. 724451 + 0. 000670* 

0. 736223-0. 007521* 

0. 439865-0. 003040*' 

0. 241064+ 0. 009891* 

0. 30 

1. 728509— 0. 003131* 

0. 735275-0. 012135* 

0. 438213-0. 003732* 

0. 239596+0. 018998* 

0.40 

1. 731756- 0. 009046i 

0.733421-0.017155* 

0. 436363-0. 003983*' 

0. 240053 + 0. 030184* 

0. 60 

1.733788-0. 016615i 

0. 730576-0. 022453* 

0. 434382 - 0. 003864*' 

0. 242847+ 0. 042932* 

0. 60 

1. 734444- 0. 025452* 

0. 726706-0. 027925* 

0. 432304-0. 003417* 

0. 248139 + 0. 066794* 

0.70 

1. 733669 - 0. 036272i 

0. 721815-0. 033484* 

0. 430143—0. 002693* 

0. 255966+0. 071449* 

0. 80 

1. 731489— 0. 045854* 

0. 715925-0. 039049* 

0. 427886- 0. 001712* 

0. 266293 + 0. 086616* 

0. 90 

1. 727947- 0. 057051*' 

0. 709075- 0. 044562* 

0. 425536-0. 000490* 

0. 279033 + 0. 102103* 

LOO 

1. 723123- 0. 068765* 

0. 701300-0. 049982* 

0. 423081 + 0. 0009681 

0. 294089 + 0. 117769* 

k 

On 

On 

On 

On 

0 

0.02 

0. 306807+0. 001570i 

0. 478591+0. 000403* 

0. 478591-0. 000403* 

0.306807-0. 001570* 

0. 06 

0. 306843+0. 004695* 

0. 478544+0. 001204* 

0. 478548- 0. 001204* 

0. 306861-0. 004695* 

0. 10 

0. 306916 + 0. 007782* 

0. 478430 + 0. 001998* 

0.478448- 0.001998*' 

0. 307002-0. 007782* 

0.20 

0. 307109 + 0. 015272i 

0. 477887+0. 003932* 

0. 478021-0. 003928*' 

0. 307779-0. 015276* 

0. 30 

0. 307572+ 0. 022405* 

0. 476691+0. 005807* 

0. 477265-0. 005767*' 

0. 309268-0. 022445* 

0.40 

0. 307077+0. 029152* 

0. 475503 + 0. 007562* 

0. 476441-0. 007496*' 

0. 311776-0. 029218* 

0. 60 

0. 306637+0. 035566* 

0. 473639+ 0. 009263* 

0. 475337- 0. 009114* 

0. 315183- 0. 035716* 

0. 60 

0. 305858+ 0. 041667* 

0. 471315+0. 010902* 

0. 474033-0. 010613* 

0. 319590-0. 041056* 

0.70 

0. 304718+0. 047474* 

0. 468635 + 0. 012481* 

0. 472541—0. 011985*' 

0. 325002- 0. 047970* 

0. 80 

0. 303222+ 0. 0530 lli 

0. 465307+0. 014008* 

0. 470852- 0. 013224* 

0.331415-0.053795* 

0. 90 

a 301382+ 0. 058287*' 

0.461644+ 0. 015488* 

0. 468970-0. 014325* 

0. 338800-0. 059460* 

1 . 00 

0. 299211+0. 063309*" 

0. 457560+0. 016922* 

0. 466894-0. 015278* 

0. 347131 — 0. 064963* 

k 

On ' 

Gn 

On 

O t , 

0 

0. 02 

0. 460172+ 0. 001446* 

0. 239286 - 0. 000715* 

-0. 239286— 0. 000717* 

-0. 460172+0. 000712* 

0. 06 

0. 460089 + 0. 002179* 

0.239261-0. 002140* 

-0. 239261-0. 002154* 

-0.460089+0. 002116* 

0. 10 

0. 459935+0. 003665* 

0.239193-0.003559* 

-0. 239193-0. 003595* 

-0. 459936 + 0. 003489* 

0. 20 

0. 459303 + 0. 007477* 

0. 238877- 0. 007074* 

-0. 238872- 0. 007208* 

-0. 459308+0. 006806* 

0. 30 

0. 458392+ 0. 011394* 

0. 238383 - 0. 010540* 

-0.238368- 0.010824* 

-0. 458407+0.009970* 

0. 40 

0. 467251 + 0. 015372* 

0. 237728- 0. 013954* 

— 0. 237693-0. 014422* 

-0. 457286 + 0. 013004* 

0. 50 

0. 455894+ 0. 019375* 

0. 236911-0. 017309 * 

- 0. 236849 - 0. 017992* 

-0. 465966+0. 016926* 

0. 60 

0. 454331 + 0. 023372* 

0. 235937-0. 020601* 

-0. 235838-0. 021518* 

-0. 454430+ 0. 018747* 

0. 70 

0. 452569 + 0. 027335* 

0. 234811-0. 023827* 

-0. 234665-0. 024981* 

-0. 452716+0. 021473* 

0. 80 

0. 450610+0. 031244* 

0. 233530-0. 026976* 

— 0. 233331-0. 028378* 

-0.450809 + 0.024110* 

0. 90 

0. 448453 + 0. 036083* 

0. 232102- 0. 030047* 

-0. 231838- 0. 031696* 

-0. 448717+ 0. 026660* 

1. 00 

0. 446108+ 0. 038835* 

0. 230529-0. 033030* 

-0. 230195-0. 034923* 

-0. 446442 + 0. 029118*' 
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TABLE I.— VALUES OF <?,»— Continued 
(b) M=0.5 


k 

<?ie 

<7,o 

0 

(?40 

0. 02 

1. 718217+0. 0008341 

0. 736335-0. 0008491 

0. 441762-0. 0005511 

0. 245279+0. 0005651 

0. 04 

1. 718021+0.0015321 

0. 736387-0. 0017211 

0. 441645-0, 0010751 

0. 244937+0. 0012651 

0. 00 

1. 719202+0. 0020881 

0. 736390-0. 0026231 

0. 441498-or 0015711 

0. 244506+0. 0021061 

0. 08 

1. 719861+0.0025051 

0.736400-0. 0035561 

0. 441332-0. 0020371 

0. 244010 + 0. 0030881 

0. 10 

1. 720582+ 0. 0027961 

a 736408-0. 0045131 

0. 441132- 0. 0024791 

0. 243471+0. 0041961 

0. 12 

1. 721371 + 0. 0029411 

0. 736408-0. 0054921 

0. 440915-0. 0028961 

0. 242899 +0. 005447! 

0. 14 

1. 722223 + 0. 0029361 

0.736379-0. 0065021 

0. 440675-0. 0032781 

0. 242316 + 0. 0068441 

0. 16 

1. 723147+ 0. 0027951 

0. 736308-0. 0076471 

0. 440405 - 0. 003672i 

0. 241733 + 0. 0085241 

0. 18 

1. 723962+0. 0026641 

0. 736244 - 0. 0085991 

0. 440155-0. 0039691 

0. 241262+0. 0100061 

0.20 

1. 724860+0. 0021771 

0. 736099-0. 0096881 

0. 439866-0. 0042741 

0. 240763 + 0. 0117851 

0.30 

1. 728804- 0. 0013381 

0. 735006-0. 0154591 

0. 438406-0. 0054961 

0. 239377+0. 0222941 

0. 40 

1. 731036-0. 0071791 

0. 732954-0. 0217271 

0.436931-0.006275! 

0. 240072+ 0. 0351811 

0. 50 

1. 732840 - 0. 0146691 

0. 729879 - 0. 0283941 

0. 435561-0. 00870541 

0. 243313 + 0. 0497681 

k 

<?., 

<7,i 

On 

<?« 

0.02 

0. 306817+0. 0020961 

0. 478688+ 0. 0005371 

0. 478584-0. 0005371 

0. 306807-0. 0020941 

0.04 

0. 306838+0. 0041821 

0. 478561 + 0. 0010721 

0. 478559-0. 0010721 

0. 306838-0. 0041821 

0.06 

0. 306892 + 0. 0062681 

0. 478498+0. 0016061 

0. 478500- 0. 0016061 

0. 306906-0. 0002581 

0. 08 

0. 306955 -i-0. 0083201 

0. 478403 +0. 0021371 

0. 478417-0. 0021301 

0. 307021-0. 0083211 

0. 10 

0. 307026+0. 0103751 

0. 478293 + 0. 0026661 

0. 478317- 0. 0026641 

0. 307160-0. 0103781 

0. 12 

0. 307102 +0. 0124071 

0. 478152+0. 0031921 

0.478201-0. 0031911 

0. 307341-0. 0124081 

0. 14 

0. 307197+0. 0144151 

0. 477988+0. 003711! 

0. 478061-0. 0037111 

0. 307550- 0. 0144151 

0. 16 

0. 307297+0. 0166051 

0. 477771 + 0. 0042801 

0. 477885 - 0. 0042791 

0. 307843-0. 0166061 

0. 18 

0. 307389+ 0. 0183811 

0. 477573+0. 0047431 

0. 477720- 0. 0047411 

0. 308114- 0. 0183821 

0. 20 

0. 307493 + 0. 0203261 

0. 477327+ 0. 0052501 

0. 477520-0. 0052491 

0. 308450-0. 0203281 

0.30 

0. 307896 + 0. 0298041 

0. 475616 + 0. 007738i 

0. 476247- 0. 0077511 

0. 311037-0. 0297921 

0. 40 

0. 308066 + 0. 0388071 

0. 473175+0. 0101231 

0. 474503- 0. 0101711 

0. 316002-0. 0387591 

0.50 

0. 307873+ 0. 0474411 

0.469953 + 0. 0124401 

0. 472474-0. 0125261 

0. 320496- 0. 0473551 

k 

<7., 

<?,, 

<?„ 

<?« 

0.02 

0. 460170 + 0. 0009581 

0. 239281-0. 0009541 

-0. 239283-0. 0009551 

-0.460168+0.0009511 

0. 04 

0. 460136 + 0. 0019241 

0. 239276-0. 0019061 

-0. 239279-0. 0019111 

-0.460133+0.0018931 

0. 06 

0. 460062+0. 0028981 

0. 239239-0. 0028561 

-0. 239245- 0. 0028701 

-0. 460066 + 0. 0028281 

0.08 

0. 459969 + 0. 0038821 

0. 239191-0. 0038061 

-0. 239191-0. 0038311 

-0. 459969+0. 0037541 

0. 10 

0. 469857+0. 0048701 

0. 239132-0.0047541 

-0. 239133-0. 004792! 

-0. 459856+0. 0046751 

0. 12 

0. 459720+0. 0058691 

0. 239058-0. 0057001 

-0. 239051-0. 0057561 

-0. 459727+ 0. 0055871 

0. 14 

0. 469570+0. 0068711 

0.238974-0. 0000431 

-0. 238967- 0. 0067201 

-0. 459577+ 0. 0064921 

0. 16 

0. 459372+ 0. 0079811 

0. 238861-0. 0076791 

-0. 238857- 0. 0077821 

-0. 459376+ 0. 0074801 

0. 18 

0.459211 + 0.008896! 

0. 238757-0. 0085281 

-0. 238763 - 0. 0086501 

-0. 459205+0. 0082821 

0.20 

0. 459017+0. 0099151 

0. 238632-0. 0094641 

-0. 238647-0. 0096191 

-0. 459002+0. 0091691 

0. 30 

0. 457749+0. 0150861 

0. 237839 - 0. 0141571 

-0. 237842- 0. 0144701 

-0. 457746+0. 0135421 

a 40 

0. 456166+0. 0203231 

0. 236761-0. 0188151 

-0. 236783-0. 0193171 

-0. 456144+0. 0178091 

0.50 

0. 454244+ 0. 0256231 

0. 235398-0. 0234591 

-0. 235440- 0. 0241591 

-0. 454202+ 0. 0219951 
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TABLE I. — VALUES OF G rn — Concluded 


(c)Af=1.0 


k 



(?J0 

Ot o 

0. 10 

2. 220453 + 0. 024422V 

0. 938021 + 0. 048128V 

0. 684779 + 0. 057930V 

0. 166747- 0. 130480V 

0. 12 

2. 219995 -i-0. 029881V 

0. 937372 + 0. 057824V 

0. 684637+0. 070004V 

0. 157996-0. 167709V 

0. 16 

2. 218702+0. 039212V 

0. 936469 + 0. 076523V 

0. 683875 + 0. 092470V 

0. 160954-0. 208206V 

0. 198675 

2. 217279+ 0. 049746V 

0. 935663 + 0. 096413V 

0. 677372+0. 117037V 

0. 169686-0. 263196V 

0. 242825 

2. 216455+0. 059393V 

0. 948201 + 0. 103049V 

0. 674169 + 0. 149419V 

0. 162176-0. 311861V 

0. 286975 

2. 212373+ 0. 070518V 

0. 930797+0. 136024V 

0. 679618+0. 16Q138V 

0. 177212- 0. 372680V 

k 

On 

On 

Ou 

0„ 

0. 10 

0. 277575+ 0. 005300V 

0. 394836+ 0. 009153V 

a 480561 + 0. 012698V 

0. 180361-0. 027161V 

0. 12 

0. 277517+0. 006386V 

0. 394828+ 0. 010933V 

0.480658+0. 015433V 

0. 180330-0. 032752V 

0. 16 

0. 277390+ 0. 008396V 

0. 394583 + 0. 014674V 

0. 480234 + 0. 020595V 

0. 181126-0. 043665V 

0. 198675 

0. 277234+ 0. 010424V 

0. 394336 + 0. 018212V 

0. 479269 + 0. 025485V 

0. 182494-0. 064121V 

0. 242825 

0. 277012+0. 012734V 

a 393992+0. 022239V 

0. 479517 + 0. 031221V 

0. 182812- 0. 066194V 

0. 286975 

0. 276750+ 0. 015034V 

0. 393575+0. 026266V 

0. 479018+ 0. 036870V 

0. 183990- 0. 078170V 

k 

On 

On 

On 

On 

0. 10 

0. 052035+ 0. 001315V 

0. 178118+0. 004118V 

0. 342280 + 0. 008298V 

0. 227667- 0. 013731V 

0. 12 

0. 052030+0. 001571V 

0. 178064+0. 004925V 

0. 342163+0. 009993V 

0. 227743- 0. 016489V 

0. 16 

0. 051996 + 0. 002091V 

0. 177993 + 0. 006633V 

0. 342012+0. 013405V 

0. 227999-0. 022129V 

0. 198675 

0. 051964+0. 002595V 

0. 177877+0. 008233V 

0. 341668+0. 016622V 

0. 228491-0. 027460V 

a 242825 

0. 051914+0. 003170V 

0. 177719 + 0. 010056V 

0. 341484+ 0. 020322V 

0. 228883 - 0. 033648V 

0. 286975 

0. 051855+ 0. 003742V 

0. 177528+0. 011878V 

0. 341115 + 0. 023999V 

0. 229602-0. 039619V 




TABLE H.— VALUES OE (7,„W FOB. THE SAMPLE CASE OF APPENDIX A 


[if-OJ; i~122J 


p 

t* 

g,«VT=? 

On vi— v 

a, ivr> 





Grl^Vl— 

ft9‘VI=? 

1 

0 

ai 

at 

a* 

0.4 

ae 

0.0 

0,7 

as 

ae 

1.721M-) 
1.71321- 
1.68706- 
1.64263-J 
1.67800-1 
1.48110- 
L37747H 
1.330M-1 
L 03310- 
0.75063- 

I-OlOOMU 
-1002231 
-aooMW 
-aoon4J 
-10020*1 
-aooiwi 
-aooisu 
-a ooieoi 
-aoois® 
-0.0000® 

aaoTss-j 

o.ioose- 

asoisi- 

129SS6- 

a 281 86- 

a 29652- 
124*02- 

a 210*1- 

a 19461- 
a 15404- 

KL 052411 

-a 012101 
-a 021201 

HI 021)81 
-a 020641 
-a 012411 
-aai7e3j 
-aaioou 

-a 013451 
-a 002771 

a 45744- 
0.45515n 

a«820- 

a 43337- 
a 41825- 
a 3831 5- 
a»565- 
a32«3- 
a 37446- 
aiw3&- 

haoiiosi 
-aoioow 
-a oi oeu 
-amo63i 
-a oioiu 
■a 00269! 
-aocessi 
-a 0078® 
■aooeeai 

-a 004811 

0 

a 017X3* 
0. 0674S- 
0.14783- 
a 25248- 
0.3727V- 
0.4959V- 
a 50253- 
O.G6118- 
0. 60783-j 


0 n 

aOC305- 
a01205H 
a 02340-j 
0.04510- 

a 06857* 
ai07«H 

a 11800- 

aiOS57H 

l-Of 

■0,00023/ 

■a 00068/ 

-0.00 IS 2/ 
-a 00322/ 
-a 00485/ 
■a 00545/ 
-a 00784/ 

-doomi 

-aOOTOli 

0 

a 00455- 

0.01783- 
0.C0937H 
0.05706- 
0.00904- 
a 13 174- 
a 10007- 
a 17555- 

a leifli- 

(-0! 

-o.oooni 

-0.000431 
-0 000261 
-0 001921 
-D, 002361 
-0 0031® 
-0 003*61 
-0 004211 
-0 00320! 

o h 

a 00017* 
aoojra 
a 01330* 
a 04040- 
a 09320- 
a 17553- 
a 30533- 
a<3315- 
149343- 

hW 
-0 i 
*0/ 

-0.00003/ 
-0.00005/ 
-a 00013/ 

-0.00023/ 
-a 00038/ 
-0.00055/ 
-0.00054J 

0 *1 
0.00003- 
0.00048-J 
0.00338- 
0.007231 
a01054n 
0.03m- 

hOI 

-01 

-0.000041 
HX 000171 
-a 00053! 
-a 001211 

-a 00*32! 
-a 003841 

-a 00*611 
-a 009411 

0 

0.00005- 
0.00073- 
0.00353- 
0.01071- 
0.02475- 
0.0474J- 
0.07843- 
0.11343- 
0. 1)033- 

h® 

-OI 

-1000021 
-1000091 
-100029! 
HXOOOCOi 
-1001141 
-1001891 
HX 002711 
-1003161 

3 

0 

ai 

0.2 

0.1 

0,4 

0,0 

0,0 

a? 

as 

ao 

0.7*7*7-1011301 
0.733*7-0. 01 1SU 
1 73247 — 0. 01107! 
0. 70341-0. 0107® 
a«76®-aO10S8i 
o.*38*e-aoo«7« 
168000-1000041 
153669-1 009071 
144242-10097® 
0. 32U1— 0. 004631 

0.47970- 
0. 47431-1 
0.49707- 
a 46474- 

145000- 

141283- 

138159- 

134043- 

0.28002- 

a20770- 

MX 009631 
-a 004901 
-a 006711 
-aoowfli 

-a 001341 
-a 005041 
-a 004*61 
-a 004101 
-a 005601 
l-a 00*641 

a 25904—0. 01066! 
a 22086- 0. 010501 
a 25323-a 010331 
a 2270S-0. 010091 
a 31317- 0. 000071 
a 50916-0. 00015! 
a 10043- 0.006441 
a 10202- 0. 0076)1 
awjss-aooesj! 
a ioj78-a oo4ooi 

0 +M 

0.007*4 -a OOOll! 
0.02820— a 000441 
0.09331 -a 00067! 

0.10813 -a 001691 

a 160*4 -a 002451 
a 21256 — a 009261 
a 26803-a 009061 
a 28316-a 004541 
a 29034 -a 003021 

0 

a 00474- 
101898- 
104063- 

aoooeiH 
a 10321- 
113722-, 
a 16981-1 
118305- 

a 10851-1 

(-01 

-aoooMi 
-a 0002® 
-a 0006® 
-a 00066! 
-a 00129! 

-aooifi® 
-a 00204! 

-a 002241 
ha 00506! 

0 +01 
aoo2*7-aoooioi 
acco33-aooo4n 
a oaow-a 000911 
a 05421-a 001661 
ao6i64-aouS28i 
a 06866- a 003041 
a 06330-000306! 
a 06140-0 00406! 
a 06405- 0. 00373! 

0 4-0* 

000007+01 
000110-0.000021 
000970- 0.0000® 
001730-0.000271 
0 03061-0.000611 
007846 -0.00117! 
012843-0.001241 
018122-0.002781 
0 21063- 0.003231 

0 

0.00005- 
0.00075- 
0.00368- 
0.01118- 
0.03580- 
0.04043- 
0.0S174- 
0.11715- 
0, 13633H 

h® 

-01 

-aoooou 

-a 00004! 
MX 000141 
■aooc3ii 
■a oooooi 
■a ooiooi 
-a 0014*1 

HX 001071 

pppppppppo 

pf&flt 

3 

0 

ai 
as 
ai 
a 4 
as 
ae 

0.7 

as 

ao 

a 44052— 0. 00484/ 
a 45831 — a 0048 Li 
a 431fl3-a 00474/ 
0.42033 -a 00452/ 
0. 40374-0. 0044X 
a 38150-100419/ 
a 38242 -a 0C2S7f 
a 31459-0. OOMW 

asoisi-aoosea/ 

a 19209-0.00211/ 

0.47703-a00«3/ 
a 47453 — a 00561/ 

0. 457)8 — a 00572/ 
0.45505 — a 00557/ 
0.4J72Q— a005W< 
0.41311— a 00 50« 

0. 35183 — a 00467* 
0.34000-a 004171 
0.29001 *-0.00350/ 
a 20783 -a 002541 „ 

-a 2J763- 0. 010731 
-a 23934- 0. 010071 
-a 23273- 0. 010611 
-a 22962-a 010231 
-a2i770-aooe8J! 
-a 20671-0. 0022BI 

-a 1000 a-0. 00969! 

-a 162*3-0. 00769! 
-a 14265-a 006441 
-110364-100468! 

0 -HW 

0.004S8-a 000061 
0.01726-a 000121 
0.0*782-1 000421 
0.0*4*0-a 000711 

a066S7-a 001061 

a 13687 -a 001321 
a i54i6-a 0017*1 
0. 10916-a 001891 

0.18664 -a 001711 

0 +0! 
ao0476-aoooo*f 
a oi870-a 000231 
ao4t»9-aooo60! 
a 09066-a 000691 
a 10328-a 00129! 
a U758-U 0010® 
a 16682- 0.00304! 
a 18817- 0.00224! 
ai*?42-aooao9! 

0 -h® 

-aooosa-aoooul 

-a 00031-0 0004*1 
-003030-0 000231 
-003488-0001571 
-a 06143-0 00232! 
-009841-000306! 
-4). 08312-0 00371! 
-006121-0 004121 
-a 06387- 0, 002 7W 

0 +01 
000004+43 
0000*2-0.000011 
0 00340 -0.000041 
0010*4-0.00001 
0 02384-0.000261 
0046*7-0.00060! 
0 075*3 - 0. 000831 

a iwB*-o. ooiioi 

012662- 0.001361 

c +® 
0.00006+® 
0.00075— a ooaoit 
0. 003fl9-Ml 000041 

o. oui»-a ooai4i 

0. 02662 — a 000311 
0. 04040 — 0. 000*01 

a osi70-a ooioo. 

0. 11723 -a 001431 
a 18643— a 001071 

0 -H2 

-0.00002+01 
-0.00037 -a 000021 

-0.00184 -a 000081 
-0.00*67 -a 000361 
-0.0128* -a 000381 
-0.05403 -a oouu 

-0.04075 -a 00184! 
-0.06887 -a 002*41 
-0.09793 -a 003071 

4 

0 

ai 

as 

aa 

a4 

ae 

ao 

a7 

a3 

ao 

a 24180- 
124066- 

aacw- 

123072- 

123197-, 

0.S0M0- 

aiB3«- 

117272-, 

ai«i3H 

0.10643- 

f-aoissw 

■1015821 

-a oiseu 
-101*2® 

-a 012731 
-0.0120*1 
-a 011121 
•loose*! 

-a 008541 

•aooooai 

130960-HX 022401 
a 30801 -a 022221 
1303*1-10218® 

a 2#sjo-a 021371 

128372-10*06® 

a 20800-a 012401 

124795-1017021 
a22io7-aoioooi 
a 18674-a 013441 
0.1*425—0 000771 

-a 45785- 
- a 45555- 
—a 44359- 
-a 41383- 
-a 41972- 
-asvaco- 

— a 96333- 
-4 33704- 
-a 37477-1 

-aiwcs- 

MX 010241 
-0.010121 
-0.01003! 
-0.00677! 
-0.00636! 
-0.009871 

-a 006161 
-a 007311 
■a 00914! 
-a 004491 

0 H 

0.003C-J 
0.00948- 
0.02070- 
0. 03547- 
0.05235- 
0.009*- 
a 05483- 
0.00283- 
0.0653V-I 

bCtf 

-a 00014/ 

-a00064i ' 

-aoouw 

-a 00103/ 

-aocaoij 
-a oo4«tf 
-a 00483/ 

-a 00634/ 

-a oo4vii 

0 +0/ 
a 00308-0, 00023/ 
aoi2i3-aoooes/ 
aoMfis-aaivu 

a 04540-0,00528/ 
a 05702- a 00485/ 
a 08915- a 0054 5/ 
a 10832- a 00784/ 
a 11887- a 00350/ 
0.10020-0,00701/ 

0 

-a 00465- 
-aai795H 
— a03932n 
-0. 05710- 
-O.OV915- 
— a 13189n 
-a 16035- 
— a 17585H 
-a 15159- 

I-CH 

■a oooio/ 
-a 00040/ 

■0. 00038/ 
-a 00150/ 
-a. 003221 
-a 00296/ 
-a 00358/ 
-a 00303/ 
-100351/ 

0 

0.00002- 
0.00038- 
a 00157* 
0,00555- 
0.01309- 
a0250S- 
a 04147- 
a 05944- 
a 05917- 

rQi 

KK 

-0.00003/ 

-0. ooou/ 

r*0. 00033/ 
-0. 00075/ 
-0. 00144/ 
rO. 00238/ 
bo. 00343/ 
-a 00396/ 

0 44W 

0. 00003+ 0! 
o.ooo«-a 000041 

0.00239 -a 00017! 
0. 00716-a 00062! 
0.01076-a 001211 
0. 03202-a 00232! 
0. 06308-a 003841 
0. 07008-a 006601 
0.06863-a 009411 

o H 
—0.00006- 
*-0.00073- 
*-0.00)54- 
—0.01075- 
-0.02475- 
-0.04748* 
-0.07852- 
-0.11254- 
-0. 1)097-1 

b 04 

-0 i 

-a 00002 / 
-a 000034 
-a 00014 / 
-a 00055/ 
-a ooioa 
-a 00175 i 
-a00«2i 

-a 002923’ 


CO 

CO 

CO 



TABLE m— VALUES OP THE DOWNWASH FACTORS Fb+Fl FOR THE NINE CONTROL POINTS OF THE SAMPLE CASE OF APPENDIX A 


[Af-OJ; t-ia] 


p 

ft 

Oontrol point 1 

Oontrol point 2 

Control point 3 

Oontrol pNnt 4 

Oontrol point 6 

Oontrol point 8 

Oontrol point 7 

Control point 8 

Control point 9 

1 

0 

0.1 

0.3 

0.3 

0.4 

0.5 

0.0 

0.7 

0.B 

0.0 

—0,311734 
—1,32030-} 
4. 07151- 
—1.20418-} 
-0.23783- 
-0.00180- 
-a 04030- 
— a 030 78- 

-acoooy 
-a 0X512- 

-4021*1! 
-409*351 
-0.203441 
-4067*41 
-0.028841 
-0,014021 
-0.010461 
-aOCBMt 
-0.000481 
NX 005594 

rfJ«J46ddddd 

T 1 1 11 1 II 1 

-0,048101 
-0 338341 
-0.842374 
-0 *29681 
-0.082371 
-0.031741 
-0.018301 
-0.011041 
-0.00691! 
-0.00714) 

—0.240744 
—1.37220-1 
3,803m- 
-1.300244 
-0.27205- 
-aiim- 
-a 05679- 
— a (4229- 

-a 02388- 
-a 02075- 

4X073*41 
4X394884 
-LC6T82) 
4X 37450) 
4X 061881 
NX 0*824) 
4X023081 
41015641 
-1011731 
-1009141 

—0.00634 1*4 
-a 017® - 
-a 05065- 
-a 03840 - 
-a 04438 - 

—a 08321 - 
-a 21380 * 
— L 24143 - 
4. 11310 - 
— L 24093 H 

4X00*9*1 
4X008001 
-0.006481 
4X007444 
4X0009*1 
-0. 012881 
-0.029641 
4X06*181 
4X206911 
4X06997) 

-101097 
-0.09188 - 
-0.0984* - 
-0.0361* - 
-0.08641 - 
-0,10*00 - 
-0.9*177 - 
-1.30166 -) 
3.9617* - 
-1.30096 -| 

4X003884 
4X008091 
4X008011 
4X010701 
-1014484 
-101860) 
-1047744 
-1292831 
-188*631 
NX 29388) 

-1011*4 4 
-103411 - 
-a 0*868 - 
-10*804 - 
-4X06968 - 
-110709 - 
-126*49 - 
-L 27778 ■ 
188715 - 
-L 27706 4 

-0.004071 
-0.010361 
-0,0117*1 
-0,0140*1 
-0.0*0831 
4X034001 
NX 07606! 
-0.364*7! 
-L0S448! 
-1364041 

-aaMEH 
-a 05714- 

-a 21028 - 
—L 24 500-1 
4.10260- 
-1.24352-} 
-a 11485- 
-0.08*21- 
-a 04328- 

NX 008861 
NX011S0! 
NX 014921 
NX 024231 
NX 064471 
NX 207881 
-106*681 
NX 029971 
nX 0138*1 
NX 00689) 

—0. 03147-1 
-0.07318- 
-0.11514- 
-0.2SCM- 
-1.30710- 
3.06755- 
-1.304144 
-a2S300- 
— a 10*00- 

-a 05617-1 

N). 0064*1 
-40168*1 
-4021741 
-4 060041 
NX 2148*! 
NXflfiltKH 
NX2S374! 
-4048181 
-401860! 
-401400! 

-a 03377- 
—0.07775' 
-a 11032- 
-0.W1O6- 
-1.18190- 
3,65250- 
-1.K0CH 
—0.28387- 
-0.10700- 
-a 05020- 

-0.0118*1 
-4 026531 
NX 018*41 
NX 07B21I 
N). *86771 
-1.06260! 
-4*8586! 
-4 07886! 
-403400! 
NX01998I 

s 

0 

0.1 

0.2 

0>3 

0.4 

as 

ao 

a7 

0.8 

0.0 

-O.Q614CH 
-Ol 11880- 
-a 1J510- 
-0.10103- 
— a 06138- 
-aassis- 
-0.02561- 

-a oi »o- 

-a 01341- 

-a 01028- 

-aoweu 

-CL 0141 a 
ha 01417/ 
4.01B31 
-aooo7i/ 
■a 00774/ 
■0. 00636/ 
-a. 00634/ 
■aooi M/ 
■0.00»8 j 

-0.1117W 
-1.12020- 
4.07351- 
-1.254 UH 
— a3^731' 
— ao^^®- 
— aO4O30- 
-0.03078- 
-0.02085- 

-a 01612- 

ha 03131/ 
H).0 «36< 
-0.202441 
hO. 007541 
-a 01554/ 
•a 01422/ 

-a 01040/ 
-aocscui 
-a 0064& 
-a 00539/ 

-a 142274 

— L 40044- 

itnoeo- 
— L 330054 
-126977- 

— 0. 11514- 
-a 06294- 
-0.03016- 
-0.02642- 
-0. 01964* 

-ao4«« 

-a 13834/ 
-164337/ 
ho. 22986/ 
ha 05*37/ 
-a 02174/ 
-0. 01630/ 
-0. 01164/ 
-0.00691 / 
-a 007141 

-0.00607 -! 
-0.01181 - 
-0.01341 - 
-0. 01W5 - 

—0. nra* * 

— 0. 03283 - 
-0.05310 - 
-a 03835 - 
-0.11000 - 
—0.06896 - 

4X009164 
4X004*7! 
4X004881 
4X004984 
NX 005884 
4X00M84 
NX 007994 
4X00900) 
-1011104 
-1009741 

-10089414 
-101780 - 
-NX 02008 - 
-NX 09840 - 
-104438 - 
-106321 - 
-111380 - 
— L *4143 - 
4. 13310 
— L 2400* 4 

-0.002031 
-1006001 
-0. 006481 
-0. 00744) 
MX 0002*1 
(-0.011381 
4X02*611 
Ml. 0*316) 
-0.20691! 
(4X062971 

-0.01027 -Z 
-0.02188 - 
-0.02042 - 
-0.03015 - 
—0.05C41 - 

-aiawe - 

-0.28177 * 
-1.301W - 
3, 0C173 
-1.30000 H 

NX 003881 
4X 008021 
NX 00901) 
-1 010701 
-1014451 
-1018691 
MX 047741 
NX 23*631 
-168*611 
1-122288) 

-0.01163-1 

-0,03000- 

-0,03816- 

-0.05713- 

-0.00*30- 

-0.11895- 

-a 00033- 
-aofoss- 
— a 0*233- 
-a 021 55- 

(4X00*161 
4). 007101 
41.007741 
4X006961 

■acumi 

4X0U86! 
4X010*51 
4X006161 
4X000451 
1-1 008311 

-4034834 
-406734- 
-409180- 
-4*1998- 
-1. *4590-1 
4.10980- 
-XJ4383! 
-4 31486- 
-4 05*21- 
-404*98-1 

-4 00580! 
-401190! 
L 4 01492! 
L 4 09423! 
MX 064471 
-4 20756! 
[-4 05*88! 
•4 022071 
■4 01266! 
[-0.00689! 

— 403147N 
-407*18- 
-411614- 
-4 18084- 
-1.30718- 
3.96758- 
-1.30414- 
-4 28300- 
-410399- 
-4 06617- 

-4008*3! 
NX 01883! 
-4021741 
NX 060041 
N> 29463! 
N) 88209! 
XO 223741 
■104318! 
■1018691 
[-101400! 

3 

0 

ai 

ao 

as 

a4 

as 

ao 

a7 

as 

0.8 

-aoio74H 
-0.02115- 
-0.02008-1 
-0. 01823- 
-0,01583- 
-0.01333- 
— O.OLlOfi- 
-O.Oa011- 
-aOO7B^- 
— a 00638- 

Ml. 002941 
NX00MS1 
-0,008441 
-0.008181 
-0. 00-1701 
-0.004041 
-0.003631 
-0.003881 
-a 003211 
1-0.002901 

-0.081481 
-411869- 
-41*510- 
-410152- 
—0, 001*8- 
— 101315- 
— 0, 02881- 
— 101809- 
-0.01341- 
-0.01028- 

MX 009611 
-a 014141 
-1014*71 
-0.012*31 
-0,009711 
-0.007741 
-0.008361 
-100834) 
-1004861 
HLOOJS6J 

-0.211714 
—1.32029- 
4.07351- 
-1.26418-1 
-a 22782* 

-aoeiaon 

-a 04930- 
-0,03078- 
-a 02096- 
-a 01512- 

ho, 0*1311 
4X09*361 
-0. 909441 
(-0.067541 
4X 028841 
4X014991 
4X010404 
4X008081 
4X006481 
4X005*91 

-0.00347 
-a 00706 - 
-a 00753 - 
-a 00027 - 
-a 00653 - 
-a 01064 - 
-a 01204 - 
-a 01312 - 
-aQ1347 - 

-a 01290 - 

1-1001681 
4X 00314) 
-100*211 
-100*311 
-100*441 
-100*001 
-0.00*781 
4X00*881 
-0.00*881 
1-0.00*741 

-0,00567 i 
—CL 01132 - 
-0.01341 - 
-0.01056 - 
-0.02220 - 
— 0. 03281 - 
—0.05310 - 
—0.09635 
-0.11000 
-0.03300 

ha 00315/ 
-a 004371 
4X00458/ 
-a 00425/ 
-a 00555/ 
-a 00048/ 
-a 00702/ 

-a oo«a/ 
-a oiiioi 

-a 00274/ 

-a00S!41h 

-a 017 ® - 

-a 02005 - 
-a 02840 * 
-a 04428 * 
-a 08321 - 
—a 21380 - 
— L M143 - 
4. 11310 
-1.14003 - 

1-0.00203/ 

-o.ooeoc/ 

-a OOC40/ 
-a 00744/ 
-a 00023/ 
-0.01288/ 
-0. 0004/ 
-O.OMIW 
-0.206011 
l-a 082971 

-a 00621 - 

—a 01207- 
-a01333h 

-a 01411 - 
-a 01490- 
-a 01493- 
-a 01406- 

-a 01240- 
-a 01004 - 
-aoosso- 

(4X009131 
4X00425! 
4X004341 
-100*40! 
-100*431 
-100*30! 
NX 00*19! 
NX0039H 
-100960! 
NX 00*98! 

—0.013C2- 

-o! 0381 5- 
-0. 05711-| 
-0.09230- 
-ansooi 
-0.09081- 
-0.05388- 
—0. 03*82- 
' — 0.0K55- 

K). 00*46! 
-0.007101 
[-0.007741 
[-0.006961 
-4010771 
-0,013881 
-0.010381 
-400616! 
-0.00048! 
-4008*1) 

-4 0J163-1 
-4 067*4- 
-0.00189- 
—4 11993- 
—1. *464)- 
4.10960- 
— L 34*8*- 
—0. 51488- 
-0. 06321- 
-0.0**2S- 

NX 00588) 
NX 01104 
-101492! 
-1094211 
-0.08447) 
-0.20788! 
MX 083861 
-0.023971 
•0.013881 
NX 00638! 

4 

0 

0.1 

0.1 

0.1 

0.4 

o.s 

0.6 

0.7 

0.8 

0.6 

-400394- 
-400790- 
HX 00783- 
-400713- 
-400061- 
-400601- 
— 400540- 

-a oweo- 

-4004**- 

-0.00*73- 

(-4001631 
-a 003231 
-aaeioi 

-0.003071 

-aooiwi 

-4001TW 
-a 003041 
-0.00*471 
-0.00*311 
-0. 002101 

-a 01074* 
-0.02115- 

-ao2006- 

-0.01823- 

-0.01683- 

-0.01333- 

-0.0UO6- 

-0.00011- 

-0.00753- 

-O.OXW- 

|NX 003841 
-1008611 
-0.008141 
-0.008181 
4X004781 
-4X004*4! 
4X003831 
4X003864 
4X003111 
4X002904 

-a 06146- 
— a 11809- 
— a 13619- 
— a 10152- 
-a«138- 
-aaJ815- 
-a 02551- 
-a 01309- 
-a 01341- 
-0. 0102$- 

(4X 009611 
-1014141 
-1014*71 
4X012*31 
4X009711 
4X 007741 
-a0083SJ 
-1006341 
-1004861 
4X00*961 

-a warn - 

-4X00415- 

-100494 

—1004*8 

— 100488 

-100474 

-1001*9 

-0.00400 

-0.0049* 

-100478 

-0.00116) 

-0.002301 

4X002311 

4X0092*1 

4X009*84 

4X009*71 

4X009*84 

4X002*71 

41.002*81 

4X002194 

-0.00347 
— aixm# 

—a 00753 
-a 00637 
-aooss3 
-a 01064 

-a 01204 
-a 01312 
-0.01347 
-a 01280 - 

(-1001604 
4X003141 
-0.00*211 
4X00*311 
■0.00*441 
•100*001 
4X 00*761 
■0.00*881 
-0.003881 
4X 00*741 

-100887 - 
-101162 - 
-101141 

-101888 

-102*28 

-10328* 

-0.06116 

-0.03038 

-111660 

—0.03896 

(NX 00218! 
NX 004*71 
NX 004861 
NX 00498! 
NXCC565! 
4X000481 
NX 007921 
-a 009001 
NX OHIO! 
-1009741 

—0.00300- 

-aooooa 

-aoowi- 

-aoo»o- 

-a 00505- 
-a 00577- 
— a 00551- 

-a 00510 - 

-a 00474- 
-0.00428- 

a ooi40/ 
■0.00280/ 
-0.00279/ 
-a 00277/ 
•0.00274/ 
-0.00167/ 
-0.00160/ 
-0.00249/ 
-0.00237/ 
-0.002241 

-0.006H- 
-0.01257- 
-a 01333- 
-a 01411- 
-a 01490- 
-a 01493- 

-a 01400- 
-a 012 -ia 
-a 01004- 
-aOC 880 - 

KL 002131 
-a 004*6/ 
•a 00434/ 
-a 00440/ 
-a 00443/ 
-a 00416/ 
-a 00419/ 
-a 00302/ 
-a oo3ao i 
-a 00318/ 

-a 01302- 

-0.02960- 
—0. 0381 B- 
-0.05712- 
-0.09230- 
-0.11895- 
-0.09083- 
-0.05188- 
—0. 03381- 
— a 02155- 

haw3«f 

-0.00710/ 

-0.00774/ 

-a 00600 / 
-0.01077/ 
-0. 01183/ 
-0.01038/ 

-a 006101 
-a 00648/ 
-a 00531/ 
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